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The  kinematic,  dynamic  and  thermodynamic  features  of  a 
simple  and  a  complex  microburst  event  were  calculated,  and  com¬ 
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Denver,  Colorado  ,.,^-difr  ing  the  Joint  Airport  Weather  Studies 
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(JAWS)  proj  ect-^  were  judiciously  edited,  gridded  and  kinemati¬ 
cally  integrated  upward  from  the  surface,  to  produce  a  three- 
dimensional  wind  field, 
three-dimensional  wind 
"goodness  of  fit"  of  the  mass  field  to  the  momentum  field**  The 
perturbation  temperature  and  thermodynamic  variables  were  used 
to  compute  the  momentum  budgets  and  field  was  derived  from  the 
buoyancy  equation,  using  the  derived  wind  and  pressure  fields. 
The  complete  set  of  Doppler  derived  winds  and  thermodynamic 
variables  were  used  to  compute  budgets  of  the  momentum  flux  and 
the  eddy  kinetic  energy  within  the  raicroburst  domain,  to  allow 
for  a  more  definitive  interpretation  of  the  results. 

>  The  results  from  this  study  reveal  some  basic  similari¬ 
ties,  along  with  some  notable  differences  in  structure.  Both 
cases  reveal  a  downward  flux  of  mass,  momentum,  and  eddying 
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The  retrieved  pressure  field  from  the 
ield  was  used  to  compute  a  relative 


energy  within  the  microburst;  the  boundary  layer  above  0.75  km 


is  the  source  of  the  mass  and  momentum;  an  intermediate  level 


near  0.7  5  km  acts  as  a  conduit  for  the  downward  flux;  and  the 


surface  layer  below  0.75  km  that  comprises  the  divergent  flow 


enclosed  in  the  microburst.  The  microburst  also  included  a  high 


reflectivity  core  at  the  center.  The  complex  case  contains  a 


misocyclone  above  0.75  km.  The  misocyclone  will  hereafter  refer 


to  the  cyclonic  circulation  above  the  surface  microburst.  The 


role  played  by  hail  was  more  significant  to  the  microburst's 


development.  The  most  significant  result  of  this  study  was  that 


a  combination  of  physical  mechanisms,  including  evaporative 


cooling,  precipitation  drag,  downward  momentum  and  eddy  kinetic 


energy  flux  and  hail  melt,  were  responsible  for  maintaining  the 


raicroburst,  during  its  mature  stage.  The  degree  to  which  any 


particular  mechanism  plays  a  role  in  the  microburst  is  depen¬ 


dent  on  the  synoptic  scale  dynamic  and  thermodynamic  struc¬ 


tures  . 
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The  downburst  was  first  identified  by  Fuj ita  (1976),  after 
a  study  of  an  aircraft  accident  at  John  F.  Kennedy  airport  in 
New  York  City,  N.Y.  in  1975.  It  was  quantified  as  a  small  down- 
burst  of  cool  air  within  a  4  km  diameter  capable  of  producing 
damaging  winds  up  to  7  5  m/ s .  The  dry  and  wet  downburst  types 
were  first  quantified  by  Fuj ita  and  Wakimoto  (1983)  in  terms  of 
measurable  precipitation  associated  with  the  downburst.  Fuj ita 
(  1985)  found  that  83Z  of  the  JAWS  microbursts  were  of  the  dry 
variety . 

A  raicroburst  was  further  defined  by  Wilson  ejt  al.  (1984) 
as  a  divergent  horizontal  outflow  of  horizontal  diameter  4  km 
with  total  horizontal  shear  >  10  m/s.  Larger  scale  phenomena 
are  referred  to  as  macrobursts.  Evaporative  cooling  was  shown 
by  Srivastava  (1985)  as  an  important  mechanism  in  initiating 
downdrafts  in  the  subcloud  layer  for  high-based  cumulonimbus 
resembling  those  found  on  the  high  plains.  The  kinematic  struc¬ 
ture  of  the  downburst  was  examined  by  Wilson  et  al.  (1984)  and 
a  pilot  of  the  dynamics  and  thermodynamics  was  completed  by  Lin 
et  al.  (1987). 

Studies  using  dual-Doppler  radar  do  not  directly  provide  a 
Cartesian  three-dimensional  wind  field,  nor  pressure  and 

1 


temperature  perturbations.  The  wind,  perturbation  pressure  and 


perturbation  temperature  fields  are  required  when  studies  of 


the  momentum  flux  and  kinetic  energy  budgets  are  made.  There¬ 


fore,  indirect  retrieval  of  these  quantities  is  required.  A 


method  for  recovery  of  the  3-D  wind  field  was  formulated  by 


Armijo  (  1969  ),  which  is  used  in  this  study.  Pressure  and  tem¬ 


perature  are  not  directly  measured  but  are  indirectly  measured 


as  the  perturbation  quantities.  These  quantities  vary  a  few 


degrees  of  temperature  or  a  few  millibars  of  pressure  from 


their  mean  value.  These  quantities  cannot  be  accurately  meas¬ 


ured  by  other  passive  or  active  sensors. 


Gal-Chen  (1978)  proposed  a  retrieval  method  where  the  per¬ 


turbation  density,  pressure  and  temperature  deviations  are 


recovered  by  solving  a  special  form  of  the  momentum  equations 


using  a  known  3-D  Cartesian  wind  field.  This  method  has  been 


used  on  both  model  and  observed  data  by  Gal-Chen  (1978),  Hane 


et  al .  (1981),  Lin  and  Pasken  (1982),  Brandes  (1984),  Roux  £t 


al .  (1984),  Hane  and  Ray  (1985),  Lin  et  al .  (1986)  and  Lin  et 


al.  (1987) 


Since  thermodynamic  retrieval  is  indirect,  data  accuracy 


is  critical  to  an  accurate  result  (Hane  et  al.,  1981) 


verification  scheme  is  required  to  insure  the  accuracy  of  the 


retrieved  variables.  Gal-Chen  and  Hane  (1981)  suggested  momen- 


mmt. 


run  checking  and  cine  continuity  as  methods.  Momentum  checking 
does  not  require  a  sequence  of  time  periods  for  comparison,  as 
time  continuity  checking  does,  so  it  is  more  widely  used. 
Momentum  checking  was  used  by  Gal-Chen  and  Kropfli  (1984),  Hane 
and  Ray  (  1985)  Lin  et  al .  (  1986)  and  Lin  et  a_l.  (1987).  Time 
continuity  checking  was  used  by  Gal-Chen  and  Kropfli  (1984). 
The  previous  studies  determined  that  the  accuracy  of  observed 
winds  was  crucial  in  determining  the  applicability  of  thermo¬ 
dynamic  retrieval. 


1_.1_  Momentum  Budgets  for  Convective  Storms 

Strong  wind  shears  in  the  boundary  layer  during  the  JAWS 
project  have  been  noted  by  Wilson  et  al.  .(1984)  and  Lin  et  al . 
(1987).  The  spatial  and  temporal  resolution  is  at  least  meso 
/-scale,  which  would  allow  for  an  investigation  of  momentum 
transfer  and  eddy  mixing  in  the  subcloud  layer.  A  study  by  Lin 
et  al .  (1987)  found  that  raicrobursts  produce  stronger  wind 
shear  and  more  effective  mixing  and  momentum  transfer. 

The  generation  of  momentum  in  the  subcloud  layer  of  a 
microburst  producing  storm  has  been  investigated  using  dual- 
Doppler  radar  by  Condray  (  1987).  He  examined  the  momentum  flux 
equation  patterned  after  LeMone  (1983)  and  the  local  flux  quan¬ 
tities.  The  results  showed  that  a  counter-gradient  momentum 


flux  occurred  below  2  km  with  strong  downdrafts  transporting 
large  positive  values  of  momentum  downward.  Additionally,  the 
buoyancy  and  pressure  terms  were  of  equal  magnitude  and  tended 
to  oppose  each  other.  The  divergence  and  shear  terms  were  weak 
secondary  contributors.  Other  authors  have  used  rawinsonde  data 
and  aircraft  measurements  to  study  the  momentum  budgets  of 
thunderstorms.  Results  indicate  that  the  horizontal  momentum 
is  approximately  conserved  following  the  parcel  in  the  updralc. 
Additionally,  the  large  momentum  produced  in  the  low-level  sub¬ 
cloud  layer  is  carried  aloft  and  a  strong  downdraft  fed  by 
potentially  cool  air  from  mid-levels  exists  (Ogura  and  Liou, 
1980).  The  vertical  transport  of  horizontal  momentum  perpendic¬ 
ular  to  a  tropical  squall  line  was  counter  the  vertical  mean 
momentum  gradient.  LeMone  (1983)  also  noted  that  the  mesoscale 
pressure  was  predominately  hydrostatic  in  the  zone  of  active 
convection.  More  research  is  needed  to  clarify  the  effect  of 
boundary  layer  raicrobursts  on  generation  of  horizontal  momentum 


wS 


1.2  Kinetic  Energy  Budgets  for  Convective  Storms 


There  are  relatively  few  studies  in'  the  literature  that 
describe  the  effect  of  areas  of  intense  convection  on  its 
immediate  environment.  Most  mesoscale,  diagnostic  kinetic 
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energy  studies  have  examined  how  large  areas  of  convection  can 
modify  the  larger  synoptic-scale  kinetic  energy  balance.  Using 
rawinsonde  data  from  the  National  Se  t  e  Storms  Laboratory 
(NSSL),  Mclnnis  and  Rung  (1972),  Rao  and  Fishman  (1974),  Rung 
and  Tsui  (  1975  )  and  Tsui  and  Rung  (  1977)  studied  the  meso  c(- 
scale  storm  environment.  Environments  containing  intense  con¬ 
vection  were  typified  by  anomalously  large  generation  and  dis¬ 
sipation  of  kinetic  energy  at  jet  stream  level.  By  contrast, 
non-convective  environments  were  characterized  by  massive  des¬ 
truction  of  kinetic  energy  from  cross-contour  flow  and  transfer 
of  energy  from  subgrid  to  grid  scales  of  motion.  Meso  a(-scale 
energy  transformations  were  an  order  of  magnitude  greater  than 
the  same  processes  on  the  synoptic  scale  for  several  environ¬ 
ments  containing  large  areas  of  intense  convection.  A  lack  of 
high  resolution  data  has  limited  research  until  recently. 

Rinetic  energy  studies  based  on  finer  resolution  data 
became  feasible  upon  completion  of  the  Atmospheric  Variability 
Experiments  (AVE)  and  the  AVE-SESAME  (Severe  Environmental 
Storms  And  Mesoscale  Experiment)  projects.  Fuelberg  and  Scog¬ 
gins  (1978)  used  AVE-SESAME  IV  rawinsonde  data  in  computing  the 
kinetic  energy  balance  of  the  mesoscale  convective  complex 
(MCC)  described  by  Maddox  (1980).  Very  large  energy  changes 
were  contained  in  areas  of  convection  embedded  in  a  rather 
inactive  region  as  a  whole.  The  large  energy  changes  in  the 


storm  environment  consisted  of  generation  or  kinetic  energy  by 


cross-contour  flow,  horizontal  flux  divergence  near  jet  stream 
level,  upward  transport  due  to  large  scale  ascending  motion  and 
transfer  of  energy  from  resolvable  to  non- resolvable  scales  of 


motion.  They  suggested  that  convection  caused  energy  changes 
since  the  maximum  of  energy  conversion  and  transport  occurred 
near  the  times  and  location  of  the  most  convection. 


Fuelberg  and  Scoggins  (1980)  used  National  Weather  Service 
(NWS)  rawinsonde  observations  in  examining  the  synoptic-scale 
energy  processes  present  during  the  Red  River  Valley  tornado 
outbreak,  which  occurred  within  the  AVE-SESAME  project.  The 
upper-level  jet  was  the  dominant  physical  mechanism  in  generat¬ 
ing  energy  transformations  at  the  synoptic  scale,  where  the 
large  magnitudes  of  horizontal  flux  convergence  and  cross¬ 
contour  destruction  of  kinetic  energy  were  associated  with 


sub-geostrophic  flow.  Conversely,  the  convective  environment 
was  dominated  by  the  generation  of  kinetic  energy  and  positive 
dissipation.  Fuelberg  et  al.  (1980)  theorized  that  sub¬ 
synoptic  scale  data  may  resolve  the  role  of  feedback  mechanisms 
between  the  storm  environment  and  the  dynamics  of  the  synoptic 
scale . 


Sub-synoptic  scale  kinetic  energy  studies  for  the  Red 
River  Valley  tornado  outbreak  were  conducted  by  Fuelberg  and 
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Jedlovec  '1  932;  and  rue  I berg  and  Printy  (1984).  The  presence  of 
additional  rawinsonde  observations  provided  finer  spatial  reso¬ 
lution  than  was  obtained  by  Fuelberg  et^  al.  (1980).  Addi¬ 
tional  results  from  Fuelberg  and  Printy  (1984)  showed  that 
the  convective  environment  was  energetically  very  active  com¬ 
pared  to  the  entire  analyzed  domain.  Furthermore,  the  magni¬ 
tudes  of  energy  transformations  in  the  convective  environment 
were  significantly  larger  and  the  effects  were  more  discernible 
than  at  the  synoptic  scale. 


A  study  of  a  modified  meso  y-scale  kinetic  energy  budget 
for  a  severe  convective  cell  embedded  within  a  squall  line  was 
completed  by  Coover  (1986).  Results  indicate  that  the  flux 
divergence  term  is  of  great  importance  within  the  convective 
cell  and  that  the  updraft/downdraft  centers  generate  the  major¬ 
ity  of  kinetic  energy  and  are  the  most  energetically  active. 
The  terms  tended  to  be  one  order  of  magnitude  greater  than  pre¬ 
vious  meso  p-scale  study  results.  The  kinetic  energy  generation 
term  was  balanced  by  the  divergence  term  rather  than  the  dissi¬ 
pation  term  of  previous  studies  at  larger  scales.  The  subcloud 
layer  was  not  included  in  the  budget  due  to  a  lack  of  data.  The 
current  research  will  concentrate  on  the  subcloud  layer. 

A  study  of  the  total  kinetic  energy  budget  conducted  by 
Rupp  (1987)  was  conducted  on  the  1647  MDT,  14  July  1982  case. 


The  results  showed  that:  1)  the  horizontal  flux  divergence  and 
horizontal  generation  terms  act  as  a  source  of  kinetic  energy 
with  turbulent  dissipation  as  a  sink;  2)  horizontal  flux 
divergence  prevails  above  5  km  AGL  with  convergence  dominating 
below.  The  vertical  and  horizontal  flux  divergence  term  nearly 
are  in  balance  and  3)  the  buoyancy  term  balances  with  the  vert¬ 
ical  generation  term. 


Pan  (1987)  conducted  an  eddy  kinetic  energy  budget  on  the 
1647  MDT,  14  July  1982  case.  Results  from  this  study,  which 
focused  on  an  ideal  case  of  a  single  microburst  with  no  signi¬ 
ficant  secondary  circulations,  demonstrate  that  eddy  kinetic 
energy  is  transported  downward,  the  horizontal  flux  divergence 
due  to  eddy  motion  acts  as  the  main  source  of  eddy  kinetic 
energy,  the  eddy  kinetic  energy  from  the  horizontal  flux  diver¬ 
gence  due  to  mean  motion  and  frictional  dissipation  act  as 
sinks  and  the  total  buoyancy  terra  is  approximately  balanced  by 
the  vertical  generation  terra. 


2.  Statement  of  the  Problem 
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Examination  of  the  dynamics  resulting  from  wet  microbursts 
is  not  presently  possible  by  means  other  than  Doppler  radar 
observations.  The  problem  is  to  process  the  dual-Doppler  data 
from  raw  signals  to  kinematic  and  thermodynamic  variables  that 
can  be  used  to  calculate  the  raicroburst  dynamics,  which 
describe  the  physical  processes  that  maintain  the  microburst. 
The  problem  is  an  inversion  problem  where  physical  variables 
present  in  the  atmosphere  such  as  temperature,  pressure  and 
three-dimensional  wind  are  recovered  from  changes  in  the 
Doppler  signal. 


The  Doppler  radar  emits  a  signal,  which  is  attenuated  by 
raindrops,  hail  and  other  obstacles  along  the  path.  The 
interaction  causes  a  phase  shift,  which  may  be  separately  cal¬ 
culated  based  upon  the  strength  of  the  return  signal.  The  col¬ 
lected  data  are  statistical  quantities,  which  require  the 
investigator  to  examine  possible  sources  of  error  and  uncer¬ 
tainty  . 


Uncertainty  in  the  data  must  be  evaluates  with  an  eye  to 
quantizing  its  impact  on  the  eventual  calculation  of  the 
dynamic  equations  and  budgets.  Because  this  research  is  fairly 
specialized  and  techniques  for  evaluating  the  impact  of 


uncertainty  are  not  always  tractable,  a  qualitative  description 


of  the  phenomenon  and  an  estimate  or  informed  speculation  may 
be  necessary. 


to 


3.  Methodology 


Doppler  radar  collected  raw  data  from  several  sites  around 
Stapleton  Airport.  NCAR's  CP-2,  CP-3  and  CP-4,  located  at  these 
sites,  formed  a  triangle  with  base  lines  of  length  15,  18  and 
28  km  (Fig.  3.1).  The  technical  characteristics  are  outlined  in 
Table  1.  The  data  were  collected  in  2  min  time  scans  by  slowly 
increasing  the  elevation  angle  per  azimuth  sweep.  A  microburst 
was  noted  at  1846  MDT  near  the  middle  of  the  analysis  domain. 
The  shallow  volume  scans  being  considered  were  centered  on 
times  1845,  1847  and  1850  MDT  for  5  August  1982. 

The  shallow  volume  scans  with  five  layer,  0.25  km  vertical 
resolution  are  ideal  for  investigating  subcloud  boundary  layer 
phenomena.  Previous  studies  by  Elmore  et  a_l.  (1986)  and  Hughes 
(1986)  on  the  5  August  1982  and  14  July  1982  cases;  respec¬ 
tively,  revealed  single  raicroburst  events  characterized  by  a 
region  of  high  pressure  at  the  downdraft  core  surrounded  by  a 
low  pressure  ring.  The  downdraft  core  was  characterized 
(Hughes,  1987)  by  strong  sinking  motion  and  warm  (positive) 
perturbation  temperature.  Induced  circulations  include  wind 
rotors  in  the  boundary  layer  on  the  upwind  side.  No  well 
defined  horizontal  cyclonic  circulations  were  discerned  or 
intensively  investigated.  A  possible  role  of  terrain  effects  on 
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3.1  Geometric  relationship  between  the  Doppler  radars 
and  the  domain  of  interest. 


Fig  • 


adar  I)i 


the  microburst  flow  has  not  been  examined  by  other  investiga¬ 


tors.  Also,  temporal  evolution  of  the  thermodynamic  variables, 
eddy  kinetic  energy,  momentum  budget  and  vorticity  budget  may 
be  examined.  Better  understanding  of  the  kinematic,  dynamic  and 
thermodynamic  structure  of  a  downburst  in  the  subcloud  layer 
over  several  minutes  is  a  goal  of  this  study. 


31.1.  The  Synoptic  Situation 

The  synoptic  situation  at  0600  MDT  on  5  August  1982 
featured  low  pressure  systems  in  southwest  Kansas  and  southeast 
Montana  (Fig.  3.2).  A  quasi-stationary  front  connected  the  low 
pressure  centers  and  bisected  northeast  Colorado.  The  front 
moved  slowly  northeast  during  the  day.  The  surface  forcing  and 
upper-level  dynamics  were  very  weak  on  this  day  with  the  jet 
stream  located  from  Oregon  through  Idaho  and  finally  exiting 
the  United  States  from  northeast  Montana. 

The  Denver  sounding  for  0000  GMT  (Greenwich  Mean  Time)  on 
5  August  1982  (Fig.  3.3)  features  a  well-mixed  boundary  layer 
extending  to  1.5  km  above  ground  level  (AGL) .  The  cumulus  con¬ 
vection  level  is  650  mb,  and  note  that  the  cumulus  temperature 
of  (30  °  C)  has  been  reached.  Also  note  the  unusual  depth  of 
the  moist  layer  (4.0  -  7.0  km)  for  a  high  plains  storm  and  the 
general  lack  of  mid-level  dry  air.  The  height  of  the  wetbulb  (0 


C)  temperature  is  2.8  km  AGL,  within  the  range  for  formation 
of  a  hail  producing  storm.  By  utilizing  the  Upper  Level  Insta¬ 
bility  (UI)  and  700  mb  dewpoint  depression  (TdD),  four  condi¬ 
tions  are  predicted  to  occur,  which  reflect  the  likelihood  of 
downburst  occurence.  With  a  700  mb  (TdD)  of  (7  °  C)  and  a  UI 
value  of  4.0,  the  predicted  region,  Area  1,  corresponds  to  an 
area  too  moist  for  downbursts,  but  thunderstorms  may  still 
develoDj Caracena  et  al.  (1983).  This  result  suggests  that  the 
dynamics  of  this  downburst  may  be  different  than  other  studies 
since  different  mechanisms  may  be  responsible  for  this  storm. 

The  synoptic  situation,  see  Fig.  (3.4),  for  the  1649  MDT, 
14  July  1982  case  is  considerably  different  than  the  5  August 
1982  case.  The  synoptic  scale  trigger  was  a  rapidly  moving  cold 
front  that  passed  through  the  domain  of  interest  at  1800  MDT. 

The  Denver  sounding,  released  at  1800  MDT,  14  July  1982  is 
graphically  presented  in  Fig.  (3.5).  The  dry  layer  extends  to  a 
height  comparable  to  the  5  August  case,  but  the  surface  mixed 
layer  has  an  average  mixing  ratio  of  4.5  g/Kg  for  the  1649  MDT, 
14  July  case,  while  a  value  of  6.5  g/Kg  occurred  in  the  5 
August  case.  The  depth  of  the  moist  layer  is  3  km  for  the  5 
August  case,  while  it  is  approximately  1.5  km  thick  in  the  14 
July  case.  Both  soundings  were  conditionally  unstable  with  lit¬ 
tle  lift  required  or  increase  of  temperature  necessary  to 


induce  convection.  The  TdD  was  (18 


o 


C)  and  the  UI  was  3.0. 


Therefore,  a  microburst  event  is  predicted  in  this  case  with 
wind  gusts  greater  than  40  knots.  Of  additional  interest  is  the 
difference  in  the  wind  profile.  The  5  August  case  features 
winds  backing  sharply  with  height  in  the  boundary  layer,  while 
the  14  July  case  has  little  directional  shear  in  the  boundary 
layer . 

3.2  Data  Reduction  and  Analysis 

Raw  Doppler  data  are  commonly  archived  for  the  researcher 
in  unified  format  (UF),  Barnes  (1973).  The  data  must  be  prop¬ 
erly  processed  to  minimize  error  and  extract  the  real  physical 
information.  Lin  and  Pasken  (1982)  outlined  the  various  stages 
in  processing  Doppler  data:  1)  Desired  data  are  extracted  and 
archived  according  to  time  and  altitude.  2)  Datum  are  con¬ 
verted  from  a  spherical  coordinate  system  to  an  orthogonal  sys¬ 
tem.  3)  Values  exceeding  pre-set  tolerances  are  eliminated  and 
"folded"  radial  velocities  are  corrected.  4)  A  spatial 
scattering  of  data  contained  in  a  given  space  is  interpolated 
to  a  two-dimensional  grid  array.  5)  The  three-dimensional  wind 
field  is  computed  from  the  gridded  fields  of  reflectivity  and 
two-dimensional  radial  velocity.  6)  Thermodynamic  retrieval  is 
performed  and  finally,  7)  These  data  are  displayed. 


The  data  consist  of  azimuth  sweeps  subtending  approxi¬ 
mately  90  deg.  of  arc.  Data  falling  outside  the  area  affecting 
the  domain  of  interest  is  deleted  by  extracting  the  proper 
sequence  of  rays.  Next,  the  data  were  converted  from  spherical 
coordinates  relative  to  the  source  radar  to  orthogonal  coordi¬ 
nates  relative  to  one  radar  (CP-2). 

The  datum  are  advected  to  a  common  time  near  the  midpoint 
of  the  data  collection  cycle  using  storm's  horizontal  motion. 
Editing  of  data  is  accomplished  using  velocity  thresholds 
predicted  by  Doppler  wave  propagation  and  a  lower  limit  on  data 
resolution  as  a  function  of  the  returned  signal.  The  resultant 
datum  were  sliced  into  slabs  centered  on  a  vertical  level  of 
interest.  The  slabs  are  then  horizontally  interpolated  to  the 
grid  points  by  use  of  the  Barnes  (1973)  scheme.  The  data  points 
are  then  vertically  interpolated  to  the  specified  horizontal 
plane . 

3.2.2  Conversion  of  Radial  Velocity  to  Cartesian  Components 
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sian  components  was  proposed  by  Armijo  (1969)  and  has  been  used 


by  Brandes  (  1977,  1978),  Lin  and  Pasken  (1982),  Lin  et  al . 


(1986)  and  Lin  et  (1987)  in  their  studies.  The  scheme  is 


governed  by  the  following  equations: 


y  _  (x~Xi)  +v(y-y,)  +  (w  j-Vt 


Eq.  (3.1) 
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Eq •  (3.2) 


where  R^  is  the  radial  distance  from  radar  1  to  target,  R  is 


the  radial  distance  from  2  to  target,  Vrl  is  the  radial  velo¬ 


city  measured  by  radar  1,  Vr2  is  the  racial  velocity  measured 


by  radar  2,  (u,v,w)  are  the  components  of  the  true  wind  in  the 


x,y,z  directions,  in  order,  (x,y,z)  is  the  actual  location  of 


the  target  with  respect  to  (w.r.t.)  the  origin,  (xl.yl.zl)  is 


the  location  of  radar  1  w.r.t.  the  origin,  (x2,y2,z2)  is  loca¬ 


tion  of  radar  2  w.r.t.  origin  and  V^_  is  the  mean  terminal  velo¬ 


city  of  hydrometeors  within  the  target  volume. 


The  terminal  fall  speed  is  estimated  using  a  formula  for  a 
high  plains  thunderstorm,  Martner  (1975),  with  a  density-height 
correction  of  Foote  and  du  Toit  (1969),  i.e., 
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Eq.  (3.3) 


where  p^(z)  is  the  environmental  air  density,  |o 


is  the  sur¬ 


face  air  density.  2  is  equivalent  radar  reflectivity  factor. 

e 

Finally,  the  anelastic  continuity  equation  is: 

SCMslw)  _  =0-  <3-‘> 

where  p^(z)  is  the  air  density  in  an  adiabatic-hydrostatic 
atmosphere . 

The  equations  in  this  section  comprise  the  equations  which 
can  be  solved  for  the  three-dimensional  wind  field  if  the  boun¬ 
dary  value  of  w  is  known.  The  continuity  equation  may  be  solved 
for  the  vertical  velocity  at  the  upper  level  and  subsequently 
integrated  upward.  The  resultant  equation  consists  of: 

WCup)  _  SaJJ-lKD  W(|Wr)  -  Eq.  (3.5) 

Vup)  (2)  Z  n 

Since  the  depth  of  the  analysis  domain  is  only  1.25  km,  the 
accumulation  of  error  is  negligible  (Lin  et  al.  ,  1986). 


2-3  Thermodynamic  Retrieval  of  Pressure  and  Temperature 

This  study  will  use  a  moving  coordinate  system.  Studies  by 
Chong  et  al.  (1980),  Hane  et  al.  (1981)  and  Gal-Chen  (1982) 
demonstrated  that  it  is  necessary  to  use  a  storm-relative  coor¬ 
dinate  system.  The  momentum  equations  used  by  Gal-Chen  (1978) 


Eq.  (3.6) 


to  compute  the  perturbation  pressure  are  as  follows: 
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where  F,  G  and  H  are  defined  as: 


Eq.  (3.7) 
Eq.  (3-8) 
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Eq. 

(3.9) 

Eq. 

(3.10) 

F,v-vw  -  -5  U 

Eq. 

(3.11) 

The  variables  include  p'  as  the  perturbation  density,  p'  is  the 
perturbation  pressure,  is  the  hydrostatic  atmosphere,  which 

is  roughly  equivalent  to  the  environmental  density  (Brandes, 
1984;  Lin  et  al .  ,  1986)  and  ^  ^ and  are  turbulent 

momentum  fluxes. 

The  turbulent  momentum  flux  terms  ,  are  computed  in  a 
manner  described  in  Appendix  A. 

From  Smagorinsky  (1963),  the  eddy  viscosity  coefficient  is 

*ui  ^ 
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Eq.  (3.12) 


Finally,  the  precipitation  drag  term,  from  Eq .  (3.11) 

is  estimated  by  relating  liquid  water  content  to  reflectivity 


as  derived  by  Kessler  (1969),  i.e., 

f  Z  (  m  ~)  '/b 
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Eq.  (3.13) 


where  b  -  1.75  and  a  -  1.73  x  10 


Once  the  values  for  precipitation  drag  and  turbulent 
momentum  fluxes  are  calculated,  Eqs.  (3.9  -  3.11)  can  be  solved 
for  the  horizontal  perturbation  pressure.  This  involves  solving 
the  following  Poisson  partial  differential  equation  resulting 
from  differentiating  Eq.  (3.6)  with  respect  to  x  and  Eq.  (3.7) 
with  respect  to  y: 

+  tZ  =  ^  Eq.  (3 


*x2 

}x 

Using  sequential  relaxation  and  the  following  Neumann 
boundary  conditions  yields: 

&  =  at 

*£'  =  Or  \  <*+  V  '  1  b 

ay 

A  unique  solution  to  Eq.  (3.14)  exists  if  the  horizontal  area 
mean  <P’>  of  P1  is  subtracted  from  P'. 


Perturbation  density  is  recovered  using  the  following 


equat ion : 

(f‘-<f>,»9  =  -(h-<h>)  ei.  (3.is: 

3  2. 

where  the  quantities  of  the  form  <>  represent  areal  means  which 
are  height-dependent  only.  The  perturbation  density,  (3'd>  and 
perturbation  pressure  are  used  to  compute  the  deviation  pertur¬ 
bation  virtual  temperature  from  the  perturbation  equation  of 
state : 


3.U  Verification  of  the  Retrieved  Thermodynamic  Data 

The  Gal-Chen  (1978)  method  is  indirect,  thus  verification 
will  be  necessary.  By  use  of  momentum  checking  and  continuity, 
Gal-Chen  and  Hane  (1981)  suggested  a  method  for  calculating  a 
correlation  factor  for  each  analysis  level.  A  non-dimensional 
factor,  Er,  is  calculated  for  the  relative  error  as  follows: 

p  -  jfi~CT£~F)  +  (fy  ~Cr)  jcfcdy  Eq.  (3.17) 

r  S$(F2  +  62)dxdy 

Gal-Chen  and  Hane  (1981)  predicted  that  values  of  Er>.5 
were  "noise",  while  Er<.5  were  scientifically  meaningful.  But, 
Hane  and  Ray  (1985)  showed  that  values  of  Er> .  5  may  still 
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contain  useful  information.  This  method  was  used  by  many 
authors.  Kropfli  (  1984),  Roux  et^  al^.  (1984),  Brandes  (1984), 
Hane  and  Ray  (  1985)  and  Lin  et  al .  (1986)  have  used  this 
method  to  verify  the  retrieved  fields.  Momentum  checking  is 
important  since  the  physical  structures  of  the  storm  should 
agree  with  fields  of  divergence/convergence  and  vertical 
motion.  Momentum  flux  reveals  the  fit  of  the  retrieved  data. 

2-5  Calculation  of  the  Turbulent  Dynamics 

Subsynoptic-scale  motion  within  convective  cells  contains 
non-hydrostatic  processes.  Synoptic-scale  kinetic  energy  budg¬ 
ets  studied  by  Kung  and  Tsui  (1975)  and  Ward  and  Smith  (1976) 
can  not  be  used  in  this  study.  This  study  considers  a  three- 
dimensional  kinetic  energy  budget  for  turbulent  processes 
occurring  in  the  surface  at  1.25  km  layer  in  a  microburst  and 
surrounding  mesoy-scale  environment. 

In  the  Doppler  data  processing,  a  three-dimensional  wind 
field  within  a  convective  cell  can  be  synthesized  from 
multiple-Doppler  data.  This  wind  field  can  be  inputed  for  cal¬ 
culating  fields  of  deviation  pressure  and  temperature  perturba¬ 
tions,  P'  and  T'd,  using  the  method  developed  by  Gal-Chen 
(  1978).  Data  resolution  for  the  JAWS  case  was  0.5  km  in  the 
horizontal  and  0.25  km  in  the  vertical.  Fine  resolution  data 
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allow  investigation  of  eddy  phenomenon  such  as  the  momentum  and 
eddy  energy  budget.  Doppler  data's  advantage  over  aircraft  data 
is  that  Doppler  data  are  nearly  simultaneously  scanned  so  the 
data  are  more  temporally  accurate. 

Values  of  the  three  Cartesian  wind  components  (u,v,w)  at 
gridpoints  can  be  used  to  calculate  the  mean  component  (u,v,w) 
and  perturbation  (u',v',w')  quantities.  The  retrieved  thermo¬ 
dynamic  perturbation  quantities  were  derived  from  their  area 
means ,  i . e . , 

P'd  -  P'  -  <P ' > ;  Tv'd  -  T'v  -  <T’v> 

not  the  environmental  mean  state.  However,  the  derived  thermo¬ 
dynamic  perturbations  can  replace  the  actual  perturbation  quan¬ 
tities  . 

3.5.1^  Calculation  of  the  Turbulent  Kinetic  Energy  Equations 

The  development  of  the  turbulent  kinetic  energy  begins 
with  the  derivation  of  the  energy  equation  for  mean  motion. 
Recall  that  the  Reynolds'  stress  tensor  is  given  by: 

XL*  =  Eq.  (3.18) 
and  the  molecular  stress  tensor  is: 

XL n  -  ~  f  <^K  + 


Eq.  (3.19) 


The  Kronecker  const  ant ,  dug  is  valued  as: 

f  i  iv  '.--k 

I  o  i** 

Additionally, 

Hik- 

deformation  tensor  =  d  Ui  .  b  Ux 

st.tic  pressure,  P  -  -1  +  ^  +  TW), 

individual  change,  D/Dt  -  and 

finally;  the  Reynolds'  equations  for  mean  motion  in  tensor  form 


is  given  by  Eq.  (3.20): 


?  Dt*  ■?  H;+  +Ti* ) 


Eq.  (3.20) 


Note  that  at  the  meso  y-scale,  the  Coriolis  term, 
2<ijkfjPUk‘  Ui“°*  B?  using  scale  analysis,  Euler's  expansion 
and  the  chain  rule^  the  equation  may  be  transformed  into: 

7  Eq.  (3.21) 

where  the  term  on  the  right  hand  side  (RHS)  of  Eq.  (3.21)  is 
the  rate  per  unit  volume  at  which  the  mean  total  apparent 
stresses  do  work  on  the  mean  motion,  and  the  tptal  rate  of 
transformation  per  unit  volume  of  mechanical  energy  of  mean 
motion  into  other  forms  of  energy,  respectively. 


WE Wf 


JO 


The  eddy  form  of  the  kinetic  energy  equation  results  from 
multiplying  Eq .  (3.21)  by  U  to  yield: 


The  terras  of  Eq.  (3.22)  can  be  expanded  by  dividing  the  total 
variable  into  a  mean  and  a  perturbed  quantity.  The  expanded 
terms  are  combined  to  form  a  new  equation.  The  equation  for  the 
mean  motion  kinetic  energy,  Eq.  (3.20),  is  subtracted  from  the 
new  equation,  which  results  in  the  eddy  kinetic  energy  equa- 


The  above  equation  is  the  equation  for  the  mean  rate  of  change 
in  the  mechanical  energy  of  the  eddying  motion,  per  unit 
volume.  The  physical  meaning  of  each  of  the  terms  is  given  as 
follows:  A  -  the  rate  of  change  of  mean  turbulent  kinetic 
energy,  per  unit  volume,  B  -  the  work  rate  of  the  geopotential 
field  on  the  eddy  flux  of  mass  arising  from  the  degree  of 
correlation  between  the  density  and  velocity  perturbation 
(i.e.,  buoyancy  production  term),  C  -  energy  cascade  effect  of 
mechanical  energy  conversion  to  the  eddy  state  from  the  mean 
state  (i.e.,  work  done  by  the  Reynolds'  stress  terms),  D  -  the 
mean  rate  at  which  the  eddying  energy  is  changing  from  work 


done  by  dissipation,  E 


the  mean  rate  of  transformation  of 


eddying  energy  into  heat  and  internal  energy,  and  F  =  the  rate 
of  Loss  of  eddying  energy,  per  unit  volume,  by  eddying  diffu¬ 
sion. 


The  molecular  stress  factor  and  Mik  are  known  and  so  sub¬ 


stitution  into  Eq.  (3.23)  yields: 

? |(i - fu r $  - f*  -  uTH  t fj an?) 
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Eq.  (3.24) 


Several  new  terms  appear  after  this  derivation.  Their  physical 
meanings  are  in  order  as  follows:  C  *  the  mean  rate  of  work  of 
the  variable  static  pressure  gradient,  D  -  the  mean  work  rate, 
per  unit  volume,  by  viscous  stresses  on  the  eddying  motion,  E  - 
the  mean  dissipation  rate,  per  unit  volume,  of  eddy  kinetic 
energy  into  heat  (>0) 


The  three-dimensional  mean  eddy  kinetic  energy,  per  unit 
volume,  may  be  written  as: 


The  buoyancy  term  (A)  can  be  rewritten  as: 


Eq.  (3.25) 
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y. 


By  recognizing: 

?' .  -t„' 


r 


l_v_ 

fv 


And 


solving  for  jo',  with  substitution  into  Eq.  (3.25),  yield: 


-puPJ  ^ 


25), 

yield : 

Eq. 

(3.26) 

in 

brackets 

represents  the  production  or  destruction  of  eddy  kinetic  energy 
due  to  precipitation  loading. 

Assume  the  mean  flow  is  a  function  of  height  only  so  that 
the  mean  values  are  area  averages  and  work  done  by  the  Rey¬ 
nolds'  stresses  are  then  described  by: 

Ti?  (3>27) 

The  molecular  dissipation  terms  can  be  grouped  together  as 
a  single  term: 

DD  =  U',)  HO.  (3.28) 

This  term  is  calculated  using  perturbations  from  the  area  means 
and  the  frictional  stresses  are  calculated  in  the  same  manner 
as  in  Appendix  A.  Substituting  Eqs.  (3.26-3.28)  into  Eq. 


(3.24)  yields  the  final  form  of  the  eddy  kinetic  energy  budget 
equation: 


The  eddy  quantities  are  calculated  in  exactly  the  same  way  as 
will  be  used  in  the  momentum  budget  calculation.  The  Reynold's 
stresses  calculated  in  the  mean  shear/turbulent  interaction 
term  were  calculated  using  the  direct  method.  The  area  average 
of  the  product  of  perturbation  wind  components  multiplied  by 
the  environmental  density  constitutes  the  direct  method. 
The  turbulent  dynamics  of  the  subcloud  layer  may  also  be  inves¬ 
tigated  in  terms  of  a  momentum  budget  that  examines  the  role  of 
eddy  and  mean  momentum  exchange.  Results  should  prove  useful  in 
confirming  results  of  the  eddy  kinetic  energy  budget  and  pro¬ 
vide  insight  into  physical  processes  not  easily  explained  by 
energy  transformations. 

3.5.2  Calculation  of  the  Vertical  Momentum  Flux 


The  local  momentum  flux  terms  are  given  by  the  following 


I  1 
-  U  VJ 


Eq.  (3.30) 


equation 


Eq.  (3.31) 


These  equations  may  be  used  to  diagnose  where  momentum  is  .tran¬ 
sported  within  the  storm.  A  vertical  slice  showing  how  the 
atmosphere  transports  available  momentum  may  be  constructed 
along  a  cross-section  of  interest. 


The  methodology  used  is  outlined  in  the  following  steps: 

1)  Storm-relative  winds  were  acquired  by  any  gridpoint 
contacting  the  cross-section. 

2)  Nearby  data  points  of  u,v  and  w  were  linearly 
interpolated  to  the  cross-section. 

3)  A  string-average  was  taken  of  all  the  points  on  the 
cross-section  unique  to  that  level. 

4)  From  the  string-averaged  mean,  perturbations  of  u',  v' 
and  w'  were  computed. 

The  location  of  the  cross-section  is  set  to  include  the 
region  of  maximum  rising  motion  and  bisecting  the  microburst. 

3.5.3  Calculation  of  the  Momentum  Budget 

Momentum  budgets  on  the  convective  scale  have  been  con- 


ducted  by  Moncrieff  and  Miller,  1976;  Sanders  and  Emanuel, 
1977;  LeMone,  1983  .  Previous  research  has  utilized  direct  tur¬ 
bulent  measurements  by  aircraft  or  model  data.  Co eve r  (1986) 
and  Condray  (1987)  utilized  retrieved  thermodynamic  variables 
from  dual-Doppler  data.  A  similar  methodology  is  employed. 

The  momentum  budget  equations  for  the  total  flow  in  a  mov¬ 
ing  Cartesian  coordinate  system  is  formulated  as: 

^  +  u|u  +V^i  =  -ij£  +£,-f»/VFx  Iq.  (3-32) 

It  +^+vli,  =  +  <3-33) 

^  +  +  +  Eq.  (3.34) 

The  Coriolis  force  is  at  least  one  order  of  magnitude 
smaller  than  the  larger  terms  of  the  above  equations.  It  is 
frequently  two  orders  of  magnitude  smaller  at  the  mesoy-scale 
and  was  not  included  in  the  budget  calculation. 

The  turbulent  component  of  the  momentum  budget  may  be 
extracted  from  the  total  flow  by  separating  it  into  a  mean  com¬ 
ponent  and  the  turbulent  (perturbed)  flow.  The  prime  value  is 
the  perturbed  part  and  the  barred  values  are  the  mean  com¬ 
ponents  of  the  area  of  interest,  which  are  calculated  from  the 
area  mean.  It  is  only  a  function  of  height.  The  variables  of 


(u , v , w, p , q )  were  likewise  separated,  which  resulted  in  the  fol 
lowing  expansion  of  Eqs .  (3.32-3.34). 


V  sLU  4-  4- 

-&b  ^ 


-ric+Ft 


Substituting  for  bP/bz  from  the  hydrostatic 

Q  +  +u!  $  +vfe)  ^  W  ^  + 

ffM®.  Wg  +'7Jte')iil+vV^£ 

=  ~y%<Th 


+ 

Eq.  (3.35) 


Eq.  (3.36) 


Eq.  (3.37) 


equation  results  i 


Eq.  (3.38) 


where  the  gravitational  buoyancy  term  can  be  rewritten  as: 


-9  ?  =9 


Eq.  (3.39) 


Schlesinger  (1980)  termed  the  quantities  on  the  right  side  of 
Eq .  (3.50)  as  the  moist  thermal  buoyancy  and  the  pressure  buoy¬ 
ancy,  respectively.  The  pressure  buoyancy  term  is  one  order  of 
magnitude  less  than  the  other  terms  so  it  is  neglected.  The  u- 
momentum  equation  (3.35)  is  multiplied  by  w'  and  the  w-momentum 
equation  is  multiplied  by  u' .  The  resultant  equations  are: 


w1  ^  -  - ufeW  =u'w’<|i4l  -  v( - 


wfeWffi  -w'W  g  + 


Eq.  (3.40) 


w1 


u‘  W  „  -  u 


_  va' W1 


(3.41) 


a  u‘(f»  -^i) 


'b*. 


Adding  Eq .  (3.40)  and  Eq.  (3.41)  and  performing  the  bar  opera¬ 


tion  yields: 


Eq.  (3.42) 
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The  terms  may  be  expanded  to  a  final  form: 

^u^.1  _  u‘  -v(z)^-®  -73SL? 

0^  V*  0  d7_  7 

- 

- - ,  \—r~r  n  — i — ,  r— >  Eq.  (3.43) 

TOw1  _W  &  +1^  -§u^ 

^  ^=2  IV 

-fwlP' 

If  ^  7 

Eq.  (3.44) 


where  term  A  «  the  local  time  change  of  u-momentum  flux, 

terra  B  -  the  mechanical  production  of  u-momentum  flux 
from  vertical  mean  shear, 

term  C  =*  a  source  of  momentum  flux  from  the 

vertical  transport  of  u-raomentum  flux, 

term  D  -  buoyancy  production  term, 

term  E  -  production  and  destruction  of  eddy  stress 

by  velocity-pressure  interactions, 

term  F  »  a  source  of  momentum  flux  from  the  horizontal 
divergence  of  u-transport, 

term  G  =*  a  source  of  momentum  flux  from  the  horizontal 
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divergence  of  v-transport, 

term  H  «  a  source  of  momentum  flux  from  the  vertical 

divergence  of  w-transport,  and 

term  I  *  the  effect  of  precipitation  loading. 

The  v-momentum  flux  may  be  derived  in  a  similar  manner  with  the 


following  result: 

®  =  -5*  %  T  ^ 

qZ. 


STS'  - 


^(uVw*)]-y[^(w  fvw^  Eq.  (3.45) 


3.5.4  Calculation  of  the  Momentum  Budgets  from 
Dual-Doppler  Data 


Momentum  budgets  have  used  either  direct  turbulent  meas¬ 
urements  using  aircraft  or  model  generated  output.  The  pertur¬ 
bation  quantities  are  not  directly  measured  by  dual-Doppler 
radar.  They  are  indirectly  measured  by  using  geometrically  syn¬ 
thesized  wind  fields,  and  assuming  an  appropriate  divergence 
profile,  to  retrieve  thermodynamic  variables.  The  retrieved 
perturbation  quantities  are  shown  to  be  equivalent  to  the 
directly  measured  variables  as  shown: 


P'd  =  P  “<P>  I  P'  =  Pd-h<P/> 

=  "C v  j  -  Tyd  ■+■ 


p  &  f  « 


(p'd ^ 


The  terms  of  the  momentum  budget  were  calculated  using 


finite  differencing  techniques.  Centered  differencing  was  used 


uO 


ir.  the  .tested  grid  of  choice  over  the  microburst,  out  some 
one-sided  differences  were  used  at  the  edges  of  the  large 
domain  of  interest. 


4.  Error  Analysis 


Error 
real  data, 
processing 
may  not  be 
An  analysis 
results  is 
research. 


analysis  and  interpretation  are  essential  when  using 
Real  data  are  subject  to  errors  during  collection, 
and  parameterization  of  physical  processes,  which 
totally  representative  of  the  real  storm  structure, 
of  possible  error  sources  and  their  impact  on  the 
therefore  necessary  to  estimate  the  validity  of  this 


Error  sources  are  of  a  random  and  non-random  nature.  Ran¬ 
dom  errors  include  ground  clutter,  inadequate  sampling,  inade¬ 
quate  temporal  synchronizing  of  the  radars,  side  lobe  contami¬ 
nation,  and  signal  processing  uncertainty.  Random  errors  are 
non-systematic ;  therefore,  it  is  virtually  impossible  to  detect 
and  correct  errors  in  the  Unified  Format  (UF)  data  because  of 


initial  data  signal  processing  and  editing.  Non-random  errors 
are  more  likely  to  be  present  in  the  data  processing  conducted 
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during  this  research,  from  the  use  of  boundary  conditions,  fin¬ 
ite  differences  methods,  and  errors  in  calculated  fields 
derived  from  errors  or  inadequately  resolved  features  in  the 
wind  field. 
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Many  researchers,  Bohne  and  Srivastava  (1976),  Doviak  et 
al .  (1976),  Ray  et  al.  (1978,1980),  Nelson  and  Brown  (1982), 
Parsons  et.  al .  (  1983  ),  Wilson  et  al..  (1984)  and  Lin,  Wang  and 
Lin  (1986),  have  quantified  random  errors  in  the  Doppler  radar 
wind  data.  All  the  researchers  concluded  that  errors  in  the 
three-dimensional  wind  field  were  <  +3  m/s.  Doviak  et.  al . 
(1976)  estimated  errors  in  the  wind  field  of  ten  degrees  in 
direction  and  1  m/s  in  speed.  Lin  and  Pasken  (1982),  Elmore 
(1982)  and  Nelson  and  Brown  (1982)  obtained  comparable  results. 
Ray  et  al.  (1978)  found  that  targets  located  within  40  km  of 
the  radars  have  a  standard  deviation  of  horizontal  wind  uncer¬ 
tainty  of  <  .2  ra/s.  Wilson  et  al.  (1984)  assumed  a  value  of 
unity  for  the  radial  velocity  variance  from  both  radars  and 
also  the  terminal  velocity  variance  to  show  that  the  combined 
error  due  to  statistical  uncertainty  in  the  radial  velocity 
estimates  and  geometric  configuration  are  approximately  1-2  m/s 
in  the  derived  horizontal  wind. 

Random  errors  in  velocity  data  are  induced  by  variance  in 
the  probability  density  of  the  assumed  drop-size  distribution. 
This  problem  is  impossible  to  quantify,  since  the  drop-size 
distribution  is  unique  to  each  storm,  even  if  it  tends  to  fol¬ 
low  a  best-fit  curve,  and  no  direct  measurements  are  available. 


-*•2  Non- random  Errors 

Non-random  errors  may  be  classified  as  error  due  to  limit¬ 
ing  assumptions  and  from  calculation  techniques,  as  outlined  by 
Hughes  (1986). 

Limiting  assumptions  include  the  choice  of  terminal 
velocity-reflectivity  relationship,  the  adiabatic-hydrostatic 
assumption  on  density,  and  the  quasi-steady  state  assumption. 


4.2^1.  Quasi-steady  State  Assumption 

Dual-doppler  radar  were  synchronized  to  scan  the  same 
storm  volume  at  the  same  time.  This  condition  was  met  for  the  5 
August  case  and  any  error  is  believed  negligible.  The  14  July 
1982  data,  where  data  collection  was  synchronized,  still  may 
have  error  since  the  radars  were  not  scanning  the  exact  space 
per  time.  The  error  would  result  as  a  difference  of  data  reso¬ 
lution  between  the  radars,  which  is  mitigated  to  a  large  extend 
by  use  of  the  Barnes'  scheme  (1973). 

In  addition  to  resolution  differences,  error  may  result 
due  to  non-linear  growth  or  decay  of  the  storm  structure  during 
the  scanning  sequence.  The  observation  time  should  be  reduced 


to  as  small  a  value  as  possible.  A  mean  correction  time  is 
assigned  as  the  center  value  of  the  scanning  sequence.  As  Gal- 
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Chen  '1932)  suggested,  a  linear  forward  or  backward  temporal 
advection  correction  is  employed.  This  condition  requires  a 
non-accelerating  storm  motion  and  a  minimum  amount  of  temporal 
change  in  storm  development. 

The  non-accelerating  storm  motion  condition  is  generally 
satisfied  as  revealed  by  Wilson  e_t  al .  (1984).  The  5  August 
1982  storm  appears  to  drift  to  the  west  at  1-2  m/s.  This  speed 
is  of  the  same  magnitude  as  the  14  July  1982  case.  The  temporal 
change  requirement  requires  an  assumption  of  constant  or  linear 
growth/decay  of  the  storm.  While  temporal  growth  is  likely  to 
be  non-linear  in  a  microburst-producing  convective  storm,  it 
apparently  does  not  introduce  unacceptable  error  for  a  two 
minute  scan  time.  This  assumption  has  been  successfully 
employed  by  Lin  and  Hughes  (1987),  Wirsing  (1985)  and  Lin  and 
Pasken  (1982). 

4 .2.2  Errors  Associated  with  the  Wind  Field 

Errors  in  this  category  include  formulation  of  the  termi- 


The  terminal  velocity-reflectivity  relationship  of  Martner 


(1975)  is  used  since  it  is  a  solution  for  a  high-plains  thun¬ 
derstorm.  Hughes  (1986)  showed  that  errors  will  be  one  order  of 
magnitude  smaller  than  the  values  of  the  u  and  v  wind  com¬ 
ponents.  He  also  showed  that  the  error  is  larger  at  greater 
elevation  angles.  The  errors  related  to  Foote  &  du  Toit  are 
also  significant  above  2  km  AGL. 

Subjective  error  in  the  wind  field  can  not  be  quantified. 
Fortunately,  few  data  points  were  in  the  14  July  1982  data  set 
and  no  subjective  correction  was  used  in  the  5  August  1982 
data.  The  method  employed  was  to  apply  restrictive  criteria  to 
the  selection  of  replacement  values.  Inconsistent  values  that 
resulted  in  large  errors  in  thermodynamic  retrieval,  were 
replaced  with  data  consistent  with  the  contoured  values  of 
non-affected  datum  points. 


4.2.3  Errors  in  the  Derived  Thermodynamic  Fields 

Errors  in  the  pressure  field  comes  from  error  in  the  wind 
field  and  the  value  of  the  mean  pressure  as  defined  by  the  data 
field.  The  effect  of  a  10Z  error  in  the  wind  field  on  the  per¬ 
turbation  field  was  examined  by  Pasken  (1981).  He  found  that 
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*he  error  in  the  wind  field  is  directly  proportional  to  the 
error  in  the  resultant  pressure  field.  The  mean  values  of  pres¬ 
sure  per  level  were  zero.  This  satisfies  the  condition  of  <p’> 
-  0,  Gal-Chen  (1978).  Lin  et  al^  (1986)  determined  that  exclu¬ 
sion  of  the  local  time  tendency  term  did  not  reduce  the  accu¬ 
racy  of  pressure  recovery. 

The  retrieved  perturbation  temperature  field  is  more  sen¬ 
sitive  to  error  than  the  perturbation  pressure  field.  Hane  et 
al .  (1981)  found  that  this  calculation  is  not  only  subject  to 
random  error,  but  may  produce  an  average  error  of  1  deg.  by 
unresolved  features  in  the  wind  field.  Results  also  indicated 
that  the  temperature  field  was  significantly  affected  by  the 
elimination  of  the  local  time  derivative. 


4.2.4  Error  in  the  Calculation  of  the  Eddy  Kinetic 

Energy  Budget  and  the  Momentum  Budget  Equations 

Positive  sources  of  error  in  these  calculations  include 
error  in  finite  difference  calculation,  applicability  of  the 
mean  value  of  the  grid  wind  field  to  the  environmental  wind  and 
the  parameterization  of  the  frictional  eddy  stresses. 

The  finite  differencing  used  was  fourth-order  on  interior 
points  with  centered  and  one-sided  differencing  used  only  on 


the  grid  edges.  The  use  of  fourth-order  differencing  maximizes 
the  data  resolution  and  minimizes  space  truncation  error. The 
mean  value  of  the  wind  was  determined  from  the  grid.  The  grid 
derived  mean  wind  was  comparable  in  magnitude  and  direction  to 
the  sounding  wind.  Little  if  any  error  is  likely  to  be  intro¬ 
duced,  because  the  5  August  1982  grid  contains  nearly  1000  grid 
points  and  the  14  July  1982  grid  contains  400  grid  points  over 
areas  of  225  and  100  km,  respectively.  Frictional  parameteriza¬ 
tion  will  introduce  error,  but  results  from  Lin  and  Condray 
(1987)  and  Pan  (1987)  suggest  that  physically  meaningful 
results  can  still  be  obtained. 

Finally,  Waldteufel  (1976)  and  Doviak  and  Zrnic  (1984) 
suggested,  and  Hughes  (1986)  showed  that  a  102  to  152  total 
error  in  the  wind  field  is  acceptable  for  accurate  results  in 
thermodynamic  retrieval.  He  also  showed  that  errors  in  reflec¬ 
tivity  were  significant  in  areas  of  high  (>  40  dBz)  areas  of 
ref lectivity. Such  an  error  would  be  of  significance  in  the  rain 
water  loading  term  of  the  perturbation-pressure  equation,  Gal- 
Chen  (1978).  However,  the  vertical  deviation  perturbation- 
pressure  gradient  term  dominates  the  buoyancy  equation  with 
rain  water  loading  only  -  102  of  the  main  term.  Therefore,  even 
a  1002  error  in  the  rain  water  loading  term  will  result  in  a 
102  error  in  the  pressure  calculation.  The  temperature  calcula¬ 
tion  will  be  more  directly  affected.  Spahn  and  Smith  (1976) 


chat  the  discrepancy  between  the  calculated  reflectivity, 


found 

based  on  an  actual  measurement  of  hail,  and  the  reflectivity 
calculated  from  the  CP-2  radar  was  in  fairly  good  agreement 
with  a  maximum  difference  of  10  -  152.  Previous  research  by 
Pasken  (1981)  showed  that  error  in  input  data,  e.g.,  wind  data, 
is  directly  proportional  in  magnitude  and  sign  to  the  thermo¬ 
dynamically  retrieved  variable.  This  result  was  verified  by 
Hughes  (1986).  For  a  more  complete  discussion  on  scale  analysis 
arguments,  see  Appendix  &. 
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The  result  of  acceptable  data  quality  and  lack  of  growth 
in  the  error  after  repeated  calculations  is  physically  realis¬ 
tic  output. 


Thermodynamic  retrieval  provided  perturbation  pressure 


fields  are  used  to  examine  the  role  of  the  three-dimensional 
pressure  gradient  in  the  maintenance  and  evolution  of  several 
wet  microbursts.  It  is  postulated  that  the  strong  outflow  and 
downdrafts  associated  with  the  microbursts  are  directly  related 
to  the  strength  of  the  perturbation-pressure  gradient.  The 
downdraft  itself  is  hypothesized  to  be  maintained  by  precipita¬ 
tion  loading  in  the  high  humidity  environment  of  a  wet  micro¬ 
burst.  Since  vertical  velocity  is  calculated  independent  of 
reflectivity,  a  positive  correlation  between  the  two  fields 
will  provide  support  to  this  theory.  The  perturbation  pressure 
fields  are  theorized  to  maintain  the  outflow  in  a  microburst. 
Eddy  kinetic  energy  budgets  and  momentum  budgets  may  be  used  as 
tracers  to  examine  how  momentum  and  kinetic  energy  transfers 
and  conversions  effect  and  maintain  the  evolution  of  micro¬ 
bursts.  The  outflow  is  postulated  to  interact  with  other  out¬ 
flows  to  form  gust  fronts,  which  may  produce  additional  convec¬ 
tive  updrafts.  A  link  between  the  microburst  outflow,  new  con¬ 
vective  updrafts  and  the  initiation  of  microbursts  will  be 
examined.  Since  two  cases,  5  August  1982  and  14  July  1982,  are 
examined,  useful  comparisons  may  be  made. 
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I*  is  theorized  that  warm  core  microburst  events  are  main¬ 
tained  by  strong  downdraft  speeds  of  more  than  7-8  m/s  in  the 
boundary  layer  in  a  wet  microburst,  due  to  the  lack  of  evapora¬ 
tive  cooling.  By  comparing  a  microburst  containing  strong  down- 
draft  speeds,  14  July  1982  case,  with  one  containing  downdraft 
speeds  of  3-4  m/s,  5  August  1982  case,  it  is  hypothesized  that 
the  weaker  case  should  contain  a  weaker  warm  core  structure  or 
perhaps  even  a  cold  core  structure.  Additionally,  microburst 
propagation  is  postulated  to  occur  as  discrete  events  organized 
by  interactions  between  a  microburst's  perturbation-pressure 
forces  and  the  updrafts  thought  to  be  associated  with  them.  By 
examining  the  case  of  a  single  microburst,  14  July  1982,  to  one 
with  multiple  microbursts,  5  August  1982,  similarities  and 
differences  in  storm  evolution  may  provide  clues  as  to  whether 
there  is  a  common  method  of  microburst  propagation. 

The  microburst  episodes  considered  occurred  on  14  July 
1982  and  5  August  1982.  Data  were  collected  from  164615-164915 
MDT  on  14  July  1982  and  1843-1851  MDT  (Mountain  Daylight  Time) 
with  data  temporally  corrected  to  1647,1649  MDT  and  1845,1850 
MDT.  Due  to  shallow  scan  heights  of  9  deg.  or  less,  vertical 
data  resolution  was  limited  to  a  2  km  AGL  (Above  Ground  Level) 
scan.  Data  analysis  and  reduction  procedures  similar  to  those 
of  Lin  et  al_.  (  1986)  and  Hughes  (  1986)  will  be  used  to  derive 
a  three-dimensional  wind  field  and  thermodynamic  variables. 
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Budge's  of  eddy  kinetic  energy  and  momentum  fluxes  at  four 
analysis  times  will  be  calculated  based  on  the  derived  wind, 
pressure  and  temperature  fields. 

By  comparing  budgets  for  different  microburst  occurrences, 
common  structures  and  processes  may  be  identified.  Differences 
may  also  be  identified  and  the  forcing  mechanisms  may  be  found. 
The  ultimate  goal  is  to  test  a  current  hypothesis  of  microburst 
structure  and  development  to  devise  a  clearer  model  of  what 
physical  processes  are  the  major  contributors  to  the  microburst 
development  and  what  circulations  may  result. 


An  analysis  of  two  separate  boundary  layers  containing 
microbursts  was  conducted  during  this  research.  The  results 
indicated  that  the  14  July  1982,  1649  MDT  case  consisted  of  one 
large  microburst  emanating  from  a  single  parent  storm;  whereas, 
the  5  August  1982,  1845  and  1850  MDT  cases  revealed  the  pres¬ 
ence  of  multiple  microburst  events.  These  two  storms  may  be 
further  categorized  as  a  quasi-steady  state  microburst,  14  July 
1982  case,  and  an  evolving  cluster  of  microbursts,  the  5  August 
1982  case,  with  the  latter  case  containing  a  microburst  in 
genesis.  Other  differences  examined  during  this  research 
include  the  role  of  terrain  and  the  possible  consequences  of 
the  presence  of  hail.  The  results  are  addressed  in  two  parts. 
The  first  part  is  an  analysis  of  the  kinematic  and  dynamic 
features  and  the  second  section  deals  with  the  dynamics 
revealed  by  the  momentum  and  eddy  kinetic  energy  budgets.  The 
following  section  will  describe  the  boundary  layer  for  the 
relatively  simple  case  of  a  single,  quas i-stat ionary  microburst 
event,  14  July  1982  case  to  the  complex  case  of  a  multiple 
raicroburst  event,  5  August  1982,  in  varying  stages  of  develop- 
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6  .  S  imple  Case :  An  Analysis  of  the  Kinemat  ics  for  the  1 4  July 
1982  Case 


The  first  microburst  event  described  will  be  the  simple 
case.  The  lack  of  interactions  between  microbursts  and  the  dom¬ 
inance  of  the  microburst  dynamics  over  the  environmental  flow 
assures  a  clearer  picture  of  the  physical  mechanisms  at  work. 
It  also  provides  a  conceptual  framework  to  compare  and  contrast 
the  complex  case  with.  The  analysis  will  be  conducted  by  exa¬ 
mining  the  horizontal  levels  for  structural  features  and  con¬ 
tinuity  of  the  features  through  the  boundary  layer. 


6.1_.1_  Elan  View  for  the  14  July  1982 ,  1649  MPT  Case 


The  plan  view  of  the  0.25  km  horizontal  wind  and  reflec¬ 
tivity  field,  Fig.  6.1,  reveals  a  symmetric  microburst  with  the 
center  located  in  the  bow  echo,  high  reflectivity  region  of 
reflectivities  greater  than  45  dBz.  A  strong  gust  front  is 
noted  northeast  of  the  center  (dashed  line),  and  a  weaker  one 
is  located  west  of  the  microburst  center.  The  strongest  winds 
occur  between  the  microburst  center  and  the  strong  gust  front 
to  the  east  of  the  center.  The  area  of  strongest  winds  east  of 
the  center  is  coincident  with  the  location  of  Standley  Lake, 
Colorado.  At  the  east  edge  of  this  area,  an  area  of  strong  con¬ 
vergence  occurs  in  a  weak  echo  region  of  reflectivities  less 


_  %•  -  *  ■ 


than  zero.  The  lack  or  signal  return  suggests  that  environmen¬ 
tal  air  is  located  there.  The  0.25  km  view  confirms  that  a 
singular  raicroburst  event  is  in  progress.  The  vertical  dimen¬ 
sions  of  the  storm  may  be  obtained  by  examining  higher  levels. 

The  0.50  km  horizontal  wind  and  reflectivity  field,  Fig. 
6.2,  reveals  that  the  core  reflectivities  have  increased  to 
values  greater  than  50  dBz.  The  horizontal  microburst  winds 
have  a  maximum  value  at  this  level,  although  the  winds  in  the 
area  of  Standley  Lake  are  of  approximately  the  same  magnitude. 
The  implication  may  be  that  the  maximum  magnitude  of  the  hor¬ 
izontal  wind  at  0.25  km  was  due  to  the  relatively  low  value  of 
friction  above  the  water  surface.  The  vertical  velocity  field 
(see  Fig.  6.3),  while  not  at  a  maximum  value,  reveals  that  twin 
upward  velocity  centers  coincided  with  the  strong  and  weak  gust 
fronts.  A  maximum  value  of  4  m/s  occurs  at  -8.5  km  east  of  CP-2 
and  -17  km  north  of  CP-2.  This  area  coincides  with  the  conver¬ 
gence  region  on  the  strong  gust  front.  A  focusing  mechanism  for 
the  updraft  location  on  the  front  may  be  the  location  of  the 
strongest  winds  in  the  cold  air  behind  the  gust  front.  A  max¬ 
imum  value  of  shear  is  encountered  where  the  maximum  winds 
intersect  the  gust  front.  The  raicroburst  driving  the  outflowing 
winds  has  a  value  of  descending  motion  of  -2  m/s  at  .5  km.  The 
raicroburst  core  does  not  contain  a  well  defined  perturbation 
pressure  bubble,  but  a  large  value  of  low  pressure  occurs  to 
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the  east  of  the  microburst  (Fig.  6.4).  The  lack  of  a  pressure 
bubble,  which  was  found  by  Lin  and  Hughes  (1987),  suggests  that 
less  mass  is  being  transported  downward  at  1649  MDT  and  more 
mass  is  being  transported  outward  with  the  stronger  horizontal 
winds  relative  to  1647  MDT,  as  the  raicroburst  enters  the  later 
stages  of  its  life  cycle.  The  large  pressure  gradient  east  of 
the  microburst  is  consistent  with  the  acceleration  of  air  to 
the  east.  Another  weaker  low  pressure  area  to  the  west  of  the 
raicroburst  coincides  with  a  weak  gust  front.  The  microburst 
center  represents  a  warm-core  system,  as  shown  in  Fig.  6.5. 
This  finding  is  consistent  with  results  obtained  by  Lin  and 
Hughes  (1987)  for  the  1647  MDT  case  that  preceded  this  one.  The 
warmest  temperatures  coincided  with  the  strong  outflow  to  the 
east  of  the  raicroburst.  The  acceleration  to  the  east  of  the 
storm  outflow  and  the  eastward  distortion  of  the  warm  anomaly, 
may  be  due  to  the  difference  in  friction  in  air  flow  over  land 
and  water  surfaces.  The  apparent  acceleration  due  to  less  fric¬ 
tion  will  yield  a  lower  pressure  at  low  levels,  where  friction 
is  a  factor.  At  a  higher  altitude  of  0.5  km,  friction  is  less  a 
factor;  therefore,  the  apparent  acceleration  due  to  the  lack  of 
friction  is  not  very  great.  The  result  of  a  lower  pressure  due 
to  the  apparent  acceleration  from  lack  of  friction,  relative  to 
the  area  away  from  the  lake,  means  that  the  dp/dz  term  in  the 
temperature  perturbation  equation  is  maximized  with  a  resultant 
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6.4  Plan  view  of  the  perturbation  pressure  field,  for 
the  14  July  1982,  1649  MDT  case  at  0.50  km. 


maximum  of  perturbation  temperature.  The  levels  from  0.75  km  to 
1.5  km  will  be  examined  using  the  reflectivity  and  horizontal 
wind  fields. 

The  horizontal  wind  flow  at  0.75  km  reveals  a  strong  gust 
front  to  the  east  of  the  microburst  is  still  intact  (Fig.  6.6). 
Of  note  is  the  appearance  of  a  dry  tongue  or  notch,  with  left 
to  right  flow  of  low  reflectivity  air,  into  the  bow  echo  struc¬ 
ture.  The  majority  of  the  horizontal  motion  contains  a  large 
component  of  mean  motion  and  comparatively  little  eddy  motion. 
Eddy  motion  consists  of  a  diffluent  pattern  over  the  bow  echo 
region  and  an  inflow  of  air  to  the  northeast  of  the  microburst 
center.  This  pattern  prevails  at  the  1.0  km  level  (Fig.  6.7), 
with  strong  winds  now  beginning  to  appear  to  the  west  of  the 
bow  echo.  Also,  strong  inflow  winds  to  the  northeast  of  the 
microburst  have  expanded  in  area  coverage  and  magnitude.  The 
trend  is  for  the  microburst  circulation  to  develop  a  jet  to  the 
rear  of  the  bow  echo  and  for  the  strong  outflow  in  front  of  the 
storm  to  lessen  in  magnitude  with  a  corresponding  dissipation 
of  the  gust  front.  The  1.5  km  horizontal  wind  pattern  reveals 
the  almost  total  loss  of  the  strong  outflow  region  east  of  the 
echo.  See  Fig.  6.8.  A  low  level  jet,  also  found  by  Lin  and 
Hughes  (1987),  at  700  mb  is  present.  Note  the  strong  west 
northwesterly  wind  maximum  to  the  west  of  the  bow  echo  at  -12 
km  east  of  CE-2  and  -19  km  north  of  CP-2.  Strong  confluence  is 
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the  north  and  northwest  of  the  wind  maximum  and  dif- 
fluence  to  the  south  and  east  of  the  wind  maximum.  The  large 
decrease  in  magnitude  of  the  horizontal  wind  across  the  bow 
echo  suggests  that  some  mass  has  to  be  transported  either 
upward  or  downward  just  southeast  of  the  50  dBz  maximum  located 
in  the  raicroburst  head.  This  phenomenon  will  be  dealt  with  in 
greater  detail  at  the  2.0  km  level. 

The  2.0  km  level  continues  the  trend  of  a  stronger  wind 
maximum  to  the  rear  of  the  microburst  (see  Fig.  6.9),  as  does 
all  levels  above  0.75  km.  The  wind  maximum  is  better  defined 
and  stronger.  It  is  interesting  to  note  that  the  area  of  max¬ 
imum  winds  is  located  in  the  area  of  maximum  reflectivity  gra¬ 
dient.  The  gradient  becomes  increasingly  tighter  in  the  conflu¬ 
ence  region  to  the  northwest  of  the  wind  maximum,  and  a 
spreading  of  the  gradient  occurs  to  the  southeast  of  the  wind 
maximum  in  the  diffluent  region.  The  vertical  velocity  field, 
shown  in  Fig.  6.10,  contains  a  large  area  of  downward  motion 
located  where  the  horizontal  winds  are  diverging  outward  from 
the  wind  maximum  across  the  bow  echo  in  the  high  reflectivity 
region.  The  reduction  in  magnitude  of  the  horizontal  wind 
across  the  bow  echo  at  the  microburst  location,  with  a 
corresponding  increase  in  downward  vertical  motion,  suggests 
some  of  the  dry,  high  speed,  environmental  air  injected  into 
the  bow  echo  mixes  with  the  saturated  downdraft  and  enhanced 


c.-.e  microburst .  Lin  and  Condray  (1988)  showed  that  momentum 
flux  is  directed  downward  for  the  1647  MDT,  14  July  1982  case. 
Strong  upward  vertical  motion  occurs  where  the  gust  front  posi¬ 
tion  occurred  at  lower  levels  to  the  area  of  the  confluence  or 
air  into  the  jet.  Apparently,  some  of  the  air  feeding  into  the 
jet  comes  from  lower  levels.  A  small  value  of  upward  motion 
associated  with  the  weak  west  gust  front  and  coincident  with 
the  horizontal  wind  maximum,  is  located  at  -19  km  north  of  CP-2 
and  -12.5  km  east  of  CP-2.  The  wind  maximum  is  associated  with 
an  area  of  lower  pressure  (Fig.  6.11).  A  strong  maximum  of  per¬ 
turbation  pressure  in  the  northeast  sector  of  the  grid 
accelerates  the  flow  to  the  southwest,  consistent  with  the  hor¬ 
izontal  wind  field.  The  temperature  perturbation  field  (Fig. 
6.12)  reveals  an  area  of  neutral  to  slightly  warm  temperature 
associated  with  the  wind  maximum,  located  in  a  low  pressure 
region.  This  field  also  reveals  an  area  of  warming  associated 
with  the  updraft  at  -8  km  east  of  CP-2  and  -16.5  km  north  of 
CP-2.  An  area  of  cooling  in  the  exit  region  of  the  wind  maximum 
is  located  at  -12.5  km  east  of  CP-2  and  -19  km  north  of  CP-2. 
This  feature  may  be  the  result  of  mixing  of  the  dry  environmen¬ 
tal  air,  advected  in  by  the  wind  maximum,  with  the  saturated 
air  of  the  precipitation  shaft.  This  area  of  maximum  cooling  is 
also  coincident  with  the  maximum  wind  crossing  the  bow  echo  at 
a  large  angle.  By  integrating  the  horizontal  plan  views  of  the 


storm  into  a  vertical  volume  and  relating  features  that  main¬ 
tain  some  continuity  with  height,  a  conceptual  model  of  the 
1649  MDT,  14  July  1982  storm,  which  represents  a  simple  case  of 
a  singular  wet  microburst,  may  be  obtained. 

The  simple  wet  microburst  in  summary  consists  of  the  clas¬ 
sic  structure  of  a  strong  downdraft  with  a  velocity  of  approxi¬ 
mately  -4  m/s  in  a  region  of  reflectivities  greater  than  45 
dBz,  spreading  out  upon  reaching  an  altitude  of  less  than  0.75 
km  AGL  (above  ground  level)  with  strong  horizontal  outflow 
below  0.5  km.  The  center  of  the  microburst  features  a  warm-core 
in  the  boundary  layer  with  a  high  pressure  perturbation  at  the 
microburst  center.  The  sole  interaction  that  takes  place  is 
between  the  microburst  and  the  surrounding  environmental  flow. 
The  environmental  flow  vector  was  given  as  the  mean  wind. 

Other  features  of  note  include  the  low  pressure  ring  that 
virtually  surrounds  the  microburst  at  low  levels  and  the  strong 
updraft  located  on  the  gust  front  to  the  east  of  the  micro¬ 
burst.  There  is  an  absence  of  strong  rotation  above  0.75  km, 
but  a  low-level  jet  (LLJ)  was  noted  for  this  case  by  Lin  and 
Hughes  (1987)  at  3.5  km.  The  low  level  jet  interacts  with  the 
downward,  saturated  air  in  the  bow  echo  to  enhance  the  down- 
draft.  Also,  this  system  was  not  rapidly  changing  between  1647 
and  1649  MDT. 
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The  simple  case  consists  of  a  rapidly  descending  stream  of 
relatively  warm  air  in  the  boundary  layer  (below  0.75  km), 
embedded  within  moderate  values  of  reflectivity  (Z  <  50  dBz) 
with  a  classic  horizontal  wind  profile  of  a  roughly  symmetric, 
circular  outflow  pattern  that  was  slowly  evolving.  The  stagna¬ 
tion  point  at  the  site  of  the  microburst's  contact  with  the 
ground  is  characterized  by  a  relatively  warm  area  of  high  pres¬ 
sure.  While  the  simple  case  may  be  advantageous  for  its  ease  of 
study  and  lack  of  complicated  interactions,  its  use  in  under¬ 
standing  how  the  microbursts  interact  with  each  other  and  the 
environment  to  promote  rapid  changes  in  the  kinematic  struc¬ 
tures  at  the  microburst  scale  (meso  -))  are  thought  to  be  too 
simplistic  for  meaningful  extrapolation.  An  examination  of  a 
more  complex  case  is  required. 


6.2  Complex  Case :  An  Analysis  of  the  Kinematics  for  the 
5  August  1982  Case 

The  second  raicroburst  event  described  will  be  the  complex 
case.  The  boundary  layer  for  the  complex  case  was  examined  at 
1845  and  1850  MDT.  The  presence  of  multiple  microbursts  and 
their  interactions  at  different  stages  of  development  allows 
the  researcher  the  ability  to  study  a  number  of  different  pos¬ 
sible  interactions  occurring  at  one  area,  at  one  time.  It  pro- 


vices  a  more  realistic  microburst  event  that  may  be  compared  t; 


the  simple  case  to  isolate  events  unique  to  the  complex  case, 
which  may  be  due  to  complex  microburst  dynamics.  Two  different 
time  periods  are  examined,  so  continuity  is  an  added  analysis 


tool  present  in  this  case. 


The  analysis  will  be  conducted  by  examining  the  horizontal 
levels  for  structural  features.  Vertical  continuity  will  be 
examined  by  evaluating  two  cross  sections  taken  through  areas 
of  interest.  The  time  continuity  of  features  is  revealed  by 


examining  two  time  periods.  The  1845  MDT  case  will  be  examined 
first  at  each  horizontal  level,  followed  by  the  1850  MDT  case. 


6.2.1.  Plan  View  for  the  5  August  1982 ,  1845  MDT  Case 


The  plan  view  of  the  0.25  km  horizontal  wind  and  reflec¬ 
tivity  field  is  shown  in  Fig.  6.13.  Distances  are  in  kilome¬ 
ters  from  the  CP-2  radar.  It  reveals  that  two  microbursts  are 
present,  Ml  and  M2,  which  occur  near  or  within  the  highest 


reflectivity  regions  of  values  greater  than  50  dBz.  The  dom¬ 
inant  microburst  is  the  Ml  feature.  Its  center  is  located  at  (x 
=  -2,  y  -  -23.5).  Note  that  the  environmental  mean  flow  at  this 


level  is  from  220  deg.  at  6-8  kts.  The  other  discernible  micro- 
burst,  M2,  is  located  to  the  northwest  of  Ml  near  (-6, -19. 5), 
and  a  third  event,  M4  (not  shown),  is  located  off  the  east  edge 
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of  the  grid  almost  due  east  of  Ml .  A  gust  front  (GF)  between 


Ml  and  M3,  marked  with  the  dashed  line,  begins  at  the  south 


margin  of  the  domain  and  extends  northward  to  two  thirds  of  the 


way  across  the  grid.  At  the  south  end  of  the  gust  front  (0, 


-25.5),  a  cyclonic  circulation  (C)  is  present.  The  wind  field 


reveals  that  the  combined  flow  of  the  environmental  wind,  and 


the  outflow  from  Ml  almost  completely  overwhelms  M2.  It  is  not 


obvious  from  the  horizontal  wind  flow  that  Ml  and  M2  are 


separate  microbursts.  The  temperature  pattern  at  this  level 


(Fig.  6.14)  may  be  of  help  in  establishing  the  fact  that  two 


distinct  microbursts  are  within  the  grid.  The  Ml  and  to  a 


lesser  extent  M2  are  both  cold-core  and  separate  entities, 


which  contrasts  with  the  warm-core  microburst  noted  for  the 


1649  MDT,  14  July  1982  case.  The  M2  microburst  does  not  have  a 


well  defined  cold  core  located  at  the  area  of  strongest  down¬ 


ward  motion  and  highest  reflectivity,  because  the  combined  out¬ 


flow  from  Ml  and  the  environment  overwhelms  and  distorts  the 


storm  structure.  Another  area  of  relatively  cool  air  is  associ¬ 


ated  with  the  gust  front.  A  very  warm  value  of  +3  C,  located 


in  the  middle  of  the  grid  at  (-2 . 5 , -18 . 5 ) ,  may  be  due  in  part 


to  terrain  and  partially  due  to  the  one  sided  finite  difference 


technique  required  at  the  lowest  level.  A  net  decrease  of  200 


meters  between  the  raicroburst  location  and  the  anomaly  location 


results  in  a  channeling  and  acceleration  of  the  flow.  An  area 
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6.14  Plan  view  of  the  perturbation  temperature  field 


for  the  5  August  1982,  1845  MDT  case  at  0.25  km. 
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of  strong  winds  at  0.25  km  is  not  reflected  at  0.50  km  because 
the  terrain  has  less  effect  at  that  level.  The  relative  low 
pressure  region  near  (-5,-17)  at  0.25  km,  associated  with  the 
area  of  strong  winds,  is  not  as  pronounced  at  0.50  km  (Fig. 
6.15),  which  yields  a  smaller  low  pressure  anomaly  relative  to 
the  surrounding  area.  The‘i^/iz  term  is  therefore  a  maximum.  A 
warm  anomaly  results  in  the  area  of  strong  downslope  winds. 

An  examination  of  the  0.50  km  level  reveals  that  the 
raicroburst  positions  on  the  grid,  and  the  cyclonic  circulation 
(C)  on  the  gust  front,  are  still  present  (Fig  6.16).  The  Ml 
microburst  still  overwhelms  the  M2  microburst  no  apparent  gust 
front  between  the  two  features.  The  Ml  microburst  features  a 
coincident  value  of  -2  m/s  vertical  motion  (Fig.  6.17).  A 
value  of  -1  ra/s  vertical  motion,  associated  with  M2,  is  removed 
from  the  region  of  highest  reflectivities,  and  is  displaced 
downwind  from  the  combined  mean  flow  and  outflow  from  Ml,  as 
the  outflow  from  Ml  overwhelms  M2 ' s  divergent  structure.  The 
maximum  divergence  is  thus  shifted  to  the  downwind  side  of  M2. 
An  important  characteristic  of  the  vertical  velocity  field  is 
that  it  clearly  delineates  Ml  and  M2.  Both  Ml  and  M2  are 
located  in  an  area  of  sinking  motion.  A  ring  of  upward  motion 
virtually  surrounds  Ml  and  M2,  even  the  area  between  Ml  and  M2. 
One  may  expect  a  correlation  between  the  vertical  velocity  and 
t.he  pressure  fields.  The  perturbation  pressure  pattern  (see 
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6.18  Plan  view  of  the  perturbation  pressure  field,  for 
the  5  August  1982,  1845  MDT  case  at  0.50  km. 
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particular  microburst .  The  microburst  downdraft  does  not 
decelerate  greatly  due  to  the  effect  of  the  earth's  surface  at 
this  level,  nor  is  it  at  a  level  where  macs  originates  and  is 
collected  for  downward  transport,  since  there  is  comparatively 
little  entrainment  into  or  transport  out  of  the  downdraft. 
Rather  it  acts  as  a  conduit,  where  mass  collected  above  may 
pass  through  this  layer  to  be  distributed  by  horizontal  diver¬ 
gence  below  0.75  km.  The  cyclonic  circulation  at  the  south  end 
of  the  gust  front  near  (-1.5, -25)  remains,  although  the  circu¬ 
lation  becomes  an  open  wave.  The  cyclonic  circulation  extends 
to  about  0.75  km  of  depth  with  ever  increasing  radius  of  rota¬ 
tion.  The  circulation  slightly  tilts  to  the  west  northwest  with 
height.  Also  note  that  the  winds  tend  to  blow  parallel  to  the 
reflectivity  contours  to  the  northeast  of  M2  and  to  the 
southeast  of  the  cyclonic  circulation.  This  suggests  that  the 
strong  outflow  of  precipitation  cooled  air  from  the  high 
reflectivity  core,  which  occurred  at  lower  levels,  is  not  tak¬ 
ing  place.  It  also  suggests  that  entrainment  of  dry  environmen¬ 
tal  air  from  flow  towards  higher  reflectivity  values  is  like¬ 
wise  absent.  The  upper  levels  should  reveal  if  any  dry  entrain¬ 
ment  of  environmental  air  occurs  at  those  levels. 


At  the  1  kra  level  (Fig.  6.21),  the  most  significant  change 
in  the  wind  flow  is  the  increase  in  the  crossing  angle  of  flow 
of  the  environmental  air  from  the  low  reflectivity  regions, 


located  southeast  of  Ml  and  northwest  of  M2,  towards  the  high 
reflectivity  cores  associated  with  the  descending  precipitation 
shaft.  This  suggests  that  entrainment  of  dry  air  is  beginning 
to  occur  at  this  level.  Also  note  the  beginning  of  cyclonic 
curvature  in  the  flow  around  the  north  sides  of  M2  and  Ml.  How¬ 
ever,  the  cyclonic  curvature  associated  with  the  cyclonic  cir¬ 
culation  on  the  south  end  of  the  gust  front  is  no  longer 
present.  The  reflectivity  pattern  features  a  finger-like 
reflectivity  pocket  southwest  of  Ml  at  0  km  east  of  CP-2.  This 
feature  curls  counter-clockwise  up  the  grid  and  back  towards 
Ml.  Another  finger-like  projection  extends  down  the  grid  from 
the  extreme  northwest  corner  of  the  grid,  pointing  at  M2.  A 
higher  level  will  be  examined  to  confirm  the  existence  of  these 
features  and  establish  a  trend  in  development  of  the  features. 

Many  structural  features  of  the  1  km  level  are  better 
organized  at  1.25  km  and  more  intense.  For  instance,  the 
fingers  of  dry  air  intrusion  are  broader  than  1  km  and  the 
crossing  angle  of  the  dry  environmental  air  toward  the  high 
reflectivity  core  increases  to  nearly  90  deg.  (Fig.  6.22).  A 
general  cyclonic  turning  of  the  winds  occurs  around  micro¬ 
bursts  Ml  and  M2.  A  distinct  cyclonic  circulation  occurs  with 
Ml  at  (-2, -23. 5).  Inflow  of  dry  air  into  Ml  is  from  the 
southeast,  while  inflow  into  M2  is  from  the  north.  Since  Ml 
effectively  scavenges  and  blocks  the  southeast  environmental 
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wind,  the  M2  microburst  circulation  is  forced  to  obtain  its 


supply  of  dry  air  from  the  north.  Unlike  Fujita's  (1985) 


description  of  the  misocyclone  structure,  where  he  described 


the  misocyclone ' s  role  as  a  collector  of  hydrometeors  for  the 


microburst,  to  enhance  precipitation  loading,  I  propose  that 


the  misocyclones  act  as  pinwheels  that  funnel  in  and  effi¬ 


ciently  collect  dry  environmental  air,  which  will  enhance  eva¬ 


porative  cooling  by  ventilating  the  saturated  downdraft.  The 


air  is  mixed  with  the  downward  rushing  air  to  create  evapora¬ 


tive  cooling  at  least  on  the  outside  layer  of  the  microburst 


and  thus  enhances  negative  buoyancy.  This  theory  is  supported 


by  looking  at  the  temperature  pattern  for  1.25  km,  as  shown  is 


Fig.  6.22.  A  warm  core  (1  C)  at  the  downdraft  center  (-2.5, 


-24.5)  with  surrounding  cool  air  (-2  C)  is  noted.  Also  note 


the  warm  (3  C)  anomaly  associated  with  the  relatively  warm, 


dry,  environmental  air  entrained  into  the  system.  More  evidence 


of  the  role  played  by  the  misocyclone  will  be  discussed  in  Sec¬ 


tion  7.  M2 ' s  circulation  is  considerably  modified  by  Ml,  espe¬ 


cially  on  its  sov*,h  side  so  comparison  is  difficult,  but  cool¬ 


ing  is  noted  on  the  north  serai-circle  of  the  microburst  and  a 


warm  anomaly  is  associated  with  the  southward  flowing  entrained 


environmental  air.  Another  warm  anomaly  at  (0,-21)  is  associ¬ 


ated  with  an  area  of  low  reflectivity  and  upward  motion  (Fig. 


6.24)  that  may  be  associated  with  the  outflow  from  the  bow 
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echo.  The  bow  echo  is  primarily  a  warm  core,  but  surrounded  by 


cooler  air,  probably  due  to  entrainment  processes.  Inflow 
environmental  air  can  be  seen  feeding  into  this  feature  from 
the  southeast.  As  the  warm  environmental  air  is  being  horizon¬ 
tally  entrained,  it  maintains  a  warm  anomaly  since  it  is  only 
very  slowly  rising  and  cooling  until  it  enters  the  main 
updraft.  Additionally,  it  does  not  encounter  evaporative  cool¬ 


ing  until  it  enters  the  updraft  and  meets  the  saturated  region 
surrounding  the  updraft.  When  it  enters  the  main  updraft  it 
also  rapidly  cools  by  dry  adiabatic  ascent;  therefore,  it 


rapidly  becomes  a  cold  anomaly.  While  the  temperature  pattern 
reveals  a  good  deal  of  the  storm's  structure,  the  pressure  pat¬ 
tern  does  not.  The  pressure  pattern  at  1.25  km  (Fig.  6.25) 
reveals  that  the  highest  pressure  is  located  in  the  southeast 
and  northwest  corners  of  the  grid.  Low  pressure  centers  are 
located  near  the  centers  of  cyclonic  circulation,  but  offset  by 
2  km.  The  mass  flow  is  therefore  toward  the  cyclonic  center 
from  the  surrounding  domain  not  associated  with  the  cyclonic 
gyre.  This  is  analogous  to  water  running  down  a  drain. 


The  areas  of  downward  motion  are  weakly  warm  at  this 
level.  Microburst  M2  is  definitely  warm-cored,  while  the  main 
updraft  is  neutral  to  slightly  cool.  One  difference  is  that  the 
distance  between  Ml  and  the  main  updraft  is  only  two  kilome¬ 
ters,  while  M2  is  more  than  four  kilometers  removed  from  the 
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updraft.  At  the  relatively  coarse  grid  scale  of  0.5  km,  the 
larger  cold  anomaly  of  the  updraft  may  influence  the  weak  warm 
anomaly  associated  with  the  downdraft.  Some  of  the  ascending 
cool  air  may  be  entrained  into  the  misocyclone,  which  in  turn 
may  add  additional  negative  buoyancy  to  the  downdraft.  The  1845 
MDT  time  period  revealed  the  presence  of  a  complex  system  of 
several  microbursts  with  associated  gust  fronts  and  raisocy- 
clones.  An  examination  of  a  second  time  period  at  1850  MDT 
will  allow  for  an  estimate  of  the  consequences  of  various 
microburst  interactions  between  themselves  "and  the  environment. 
A  conceptual  model  may  be  formulated  explaining  how  the  storms 
are  evolving. 


6.2.2  Plan  View  for  the  5  August  1982 ,  1850  MDT  Case 

The  plan  view  of  the  0.25  km  level  will  be  examined  by 
looking  at  the  0.50  km  grid  spacing  case  and  one  that  uses  0.25 
km  grid  spacing.  The  0.50  km  grid  (Fig.  6.26)  shows  that  the 
wind  circulation  pattern  undergoes  some  subtle  but  significant 
changes  by  1850  MDT.  The  major  change  is  the  appearance  of  M3 
at  the  eastern  tail  end  of  the  bow  echo,  associated  with  Ml. 
The  50  dBz  contour  has  expanded  and  M3's  appearance  is  coin¬ 
cident  with  the  appearance  of  the  50  dBz  contour.  The  gust 
front  is  very  well  defined,  but  it  is  bridged  by  M3  at  its  raid- 
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Fig.  6.26  Plan  view  of  the  horizontal  wind  and  reflectivity 
field,  for  the  5  August  1982,  1850  MDT  case  at 

0.25  km. 
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die  point.  M3  has  aiso  succeeded  in  disrupting  the  inflow  into 
the  cyclonic  gyre  on  the  gust  front.  The  disruption  of  the 
inflow  into  the  cyclonic  gyre  on  the  gust  front  has  forced  the 
cyclonic  gyre  southeast  along  the  gust  front  to  the  edge  of  the 
grid.  The  movement  south  results  from  the  southward  shift  in 
the  area  on  the  gust  front  that  interacts  with  the  southeast 
environmental  flow.  The  area  of  greatest  cyclonic  shear  com¬ 
bined  with  interaction  with  warm,  buoyant  environmental  air 
occurs  at  this  point.  The  gyre  is  in  a  weaker  state  of  organi¬ 
zation.  The  M2  microburst  remains  embedded  within  the  outflow 
of  Ml  and  the  environmental  flow;  however,  the  M2  microburst 
does  possess  greater  values  of  reflectivity  at  1850  MDT  with 
values  of  55  dBz  present.  There  is  a  strong  flow  of  air  toward 
lower  dBz  values  of  reflectivity.  This  suggests  that  precipita¬ 
tion  cooled  air  is  flowing  toward  and  undercutting  the  warmer 
environmental  air.  Recall  from  the  discussion  of  the  1845  MDT 
plan  views,  at  higher  levels,  the  flow  is  toward  the  high 
reflectivity  values. 

This  level  will  now  be  examined  at  the  smaller  grid  scale 
of  0.25  km  (Fig.  6.27)  to  see  what  sensitivity  the  main 
features  may  have  to  a  higher  resolution  grid.  One  difference 
is  that  the  Ml  bow  echo  now  has  reflectivities  greater  than  60 
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Fig.  6.27  Same  as  Fig.  6.26,  except  using  a  0.25  km  grid 
spac ing . 


it  is  composed  of  several  distinct  cells;  whereas,  at  the  0.50 
km  scale  it  was  rather  amorphous.  It  is  interesting  to  note 
that  all  three  of  the  microbursts  are  embedded  within  an  area 
of  reflectivity  greater  than  55  dBz.  Miller  et  al^  (1988)  used 
CP-3  and  CP-4  radars  to  investigate  a  severe  Colorado  super¬ 
cell.  They  noted  that  locations  of  hail  reports  were  mostly 
inside  the  swath  of  the  55  dBz  contour.  The  M2  microburst 
reveals  that  the  reflectivity  pattern  appears  to  be  evolving  a 
bow  shaped  configuration.  The  general  airflow  appears  to  have  a 
stronger  northerly  component  to  the  east  of  CP-2  than  the  0.50 
km  case.  As  in  the  0.50  km  case,  the  areas  of  strongest  outflow 
are  predominately  on  the  north  side  of  the  microbursts,  where 
the  environmental  winds  and  the  microburst  winds  sum.  Also  of 
interest  is  the  confluence  zone  to  the  east  northeast  of  the  60 
dBz  maximum.  This  feature  appears  to  be  channeled  flow  caused 
by  the  interaction  of  the  bow  echo  microburst  complex  with  a 
microburst  located  to  the  east  of  the  grid.  This  feature  is 
much  better  defined  at  the  smaller  scale.  The  pressure  pattern 
shown  in  Fig.  6.28,  maintains  the  same  general  pressure  pattern 
as  for  the  0.50  km  case,  but  pressure  gradients  are  generally 
20  to  25Z  greater.  The  most  striking  change  is  the  large  pres¬ 
sure  gradients  that  occur  on  the  north  side  of  the  outflow 
areas.  These  areas  reflect  the  underlying  strong  small  scale 
wind  fields.  These  wind  fields  are  not  resolvable  at  a  0.50  km 


6.28  Plan  view  of  the  perturbation  pressure  field,  for 
the  5  August  1982,  1850  MDT  case  at  0.25  km,  for 
a  0.25  km  grid  spacing. 
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scale  cue  to  che  inherent  smoothing  of  the  Barnes  (1973) 


scheme.  The  gradient  in  the  area  north  of  M2,  between  the 


microburst  and  the  maximum  outflow  winds,  approaches  1  mb 


between  M2  and  (-6,  -17.5).  Along  the  gust  front,  a  sharp  pres¬ 


sure  gradient  and  low  pressure  mark  the  location  of  maximum 


convergence.  This  same  feature  has  a  value  of  upward  motion  of 


2  m/s  for  the  0.25  km  case  (Fig.  6.29);  whereas,  the  0.50  km 


case  had  a  value  of  1  m/ s .  Dowcward  vertical  velocity  of  -1.5 


m/s  is  associated  with  Ml  for  the  0.25  km  grid  spacing.  The 


value  of  -1.0  m/s  at  the  0.5  km  grid  scale  suggests  that  the 


updraft  is  more  sensitive  to  the  change  in  resolution,  probably 


because  the  enhanced  resolution  of  the  high  speed  outflow  near 


the  outer  edges  of  the  microbursts  produces  enhanced  conver¬ 


gence  at  the  gust  front  boundaries.  Additionally,  the  finer 


resolution  allows  for  sharper  definition  and  greater  discon¬ 


tinuity  of  the  wind  flow  near  the  gust  front,  which  will  result 


in  greater  values  of  convergence  and  upward  motion.  Other 


features  include  the  confluence  line  to  the  northeast  of  the 


bow  echo,  associated  with  upward  motion  of  0.5  m/s,  and  the 


extreme  upper  left  corner  of  the  grid  is  an  area  of  upward 


motion  associated  with  a  gust  front  emanating  from  M2. 


The  0.50  km  wind  pattern,  shown  in  Fig.  6.30,  is  essen¬ 


tially  the  same  pattern  as  occurred  at  the  0.25  km  level.  Once 


again,  the  areas  of  strongest  outflow  are  on  the  north  side  of 


Fig.  6.29  Plan  view  of  the  vertical  motion  field,  for  the  5 
August  1982,  1850  MDT  case  at  0.25  km,  for  a  0.25 
km  grid  spacing. 
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6.30  Same  as  Fig.  6.26,  except  for  0.50  km. 


■  r.e  micrcbursts ;  however,  between  the  outflow  from  Ml  and  the 
microburst  M2,  the  winds  are  much  reduced  in  magnitude.  This 
suggests  that  M2  tries  to  oppose  the  combined  Ml  and  environ¬ 
mental  flow.  The  microburst  M2  has  its  greatest  intensity  at 
this  level.  The  wind  flow,  in  general,  is  at  a  maximum  at  0.50 
km  due  to  the  lessened  effect  of  surface  friction  and  the 
better  organization  of  the  microburst.  The  vertical  velocity 
field  (see  Fig.  6.31)  features  the  bow  echo  containing  Ml  and 
M3,  having  a  maximum  value  of  -2  m/s  near  the  Ml  microburst.  An 
updraft  velocity  of  2  m/s  at  (0,-25)  is  associated  with  the 
gust  front  between  Ml  and  M3.  This  value  is  double  the  value 
for  1845  MDT.  The  downdraft  Ml  is  the  same  magnitude  as  1845 
MDT,  but  M3  is  assisting  in  making  the  bow  echo  an  area  of 
downward  motion.  A  semi-circle  of  downward  motion  surrounds  the 
updraft.  The  intensification  of  convergence  and  upward  motion 
at  the  gust  front  appears  due  to  M3's  intensification.  The  M3 
microburst  will  create  a  vigorous  outwardly  propagating  gust 
front  that  occludes  with  the  old  gust  front.  At  the  triple 
point  of  the  occlusion,  the  weakened  cyclonic  circulation  is 
found.  While  the  microburst  M2  is  not  obvious  from  the  horizon¬ 
tal  wind  pattern,  it  is  more  obvious  from  the  vertical  motion 
pattern.  Notice  how  M2  is  separated  from  Ml  and  M3  by  an  area 
of  upward  motion  between  the  two  systems.  M2  was  not  classified 
as  a  separate  system  by  storm  cataloguers  based  on  the  horizon- 


cal  wind  shear  requirements  for  a  microburst.  It  is  apparent 
from  the  vertical  velocity  field  that  two  distinct  centers  of 
downward  motion  are  present.  Also,  two  distinct  centers  of  50 
dBz  plus  reflectivities  are  noted.  The  result  suggests  that  the 
present  convention  of  classifying  microbursts  is  based  on  the 
singular,  isolated  microburst  event.  A  broader  definition  based 
on  the  total  storm  structure  to  include  the  vertical  velocity 
field  may  prove  more  accurate.  The  perturbation  pressure  pat¬ 
tern  at  0.50  (Fig.  6.32),  reveals  that  the  highest  value  of 
perturbation  pressure  is  located  on  the  upwind  side  (to  the 
southeast)  of  the  microburst.  The  pressure  pattern  reveals  that 
Ml,  M2  and  M3  are  distinct  systems.  The  maximum  pressure  gra¬ 
dient  for  this  case  exists  between  M2  and  the  low  pressure  to 
the  north  at  (-5 . 5 , -16 . 5) .  A  value  of  0.2  mb/km  is  observed.  If 
tendency  is  included,  the  value  of  the  pressure  gradient  is 
increased  to  .4  rab/kra,  not  shown.  The  inclusion  of  tendency 
increased  the  amplitude  of  the  microburst  pressure  field.  The 
temperature  structure  should  further  define  the  microburst 
structure.  The  temperature  pattern  at  0.50  km  (Fig.  6.33) 
reveals  that  the  gust  fronts  are  cool  anomalies  due  to  upward 
motion  and  adiabatic  expansion  found  in  the  updrafts.  The  down- 
drafts  are  also  cool  at  1850  MDT,  as  they  were  at  1845  MDT.  The 
14  July  1982  case  featured  a  warm-core  microburst.  The  struc¬ 
tural  difference  that  accounts  for  part  of  the  difference 
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between  the  two  cases  seems  to  be  that  the  5  August  1982  case 
features  slower  descending  downdrafts,  which  allows  the 
environmental  dry  air  and  the  descending  saturated  air  to  mix 
and  evaporat ively  cool.  An  additional  cause  for  the  discrepancy 
may  include  cooling  from  entrainment  of  the  dry  environmental 
air  above  0.75  km,  which  may  induce  additional  cooling  of  the 
updraft.  To  fully  understand  the  role  played  by  the  raisocyclone 
structure,  which  occurs  above  the  0.75  km  surface  layer  micro¬ 
burst  layer,  the  higher  levels  will  be  examined  in  turn. 

The  0.75  km  wind  and  reflectivity  pattern  (Fig.  6.34), 
shows  the  dominance  of  the  environmental  flow.  As  during  the 
1845  MDT  case,  the  flow  in  the  inflow  region  southeast  of  Ml 
and  northwest  of  M2  is  almost  parallel  to  the  reflectivity  con¬ 
tours,  except  at  the  northern  fringe  of  the  M2  inflow.  The  gust 
front  manages  to  maintain  its  structure  at  this  height; 
whereas,  at  1845  MDT,  it  had  largely  disappeared.  The  gust 
front  apparently  is  enhanced,  even  at  higher  levels,  by  the 
growth  of  the  M3  raicroburst.  Unlike  at  1845  MDT,  the  raisocy¬ 
clone  has  appeared  at  the  head  of  the  bow  echo  at  (-3, -24. 5) 
near  Ml.  The  descent  of  the  misocyclone  to  0.75  km  from  1  km 
suggests  that  the  raicroburst  system  has  intensified.  The 
greater  reflectivity  present  suggests  that  as  the  reflectivity 
core  descends,  and  the  system  spins  up,  the  raisocyclone  can 
extend  to  lower  levels  because  the  spin  up  compresses  the 
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6.35  Same  as  Fig.  6.26, 


except  for  1.0  km 


shows  that  the  misocyciones  are  now  very  well  organized.  The 
head  of  the  Ml  bow  echo  contains  a  well-defined  misocyclone  at 
(-3,-24. 5).  Even  the  M2  microburst  has  developed  a  well  defined 
misocyclone,  to  its  right  at  (-4. 5, -19).  The  increased  organi¬ 
zation  of  the  raicrobursts,  coupled  with  the  increased  value  of 
downward  motion  and  the  previously  noted  descent  and  compres¬ 
sion  of  the  misocyclone  suggests  that  the  spin  up  of  the  miso¬ 
cyclone  parallels  the  increase  in  intensity  of  the  surface 
layer  microburst.  The  microburst  systems  are  therefore  intensi¬ 
fying.  Additionally,  the  source  region  of  dry  environmental 
inflow  into  M2  now  includes  environmental  air  from  the 
northeast  of  the  misocyclone  circulation.  The  reflectivity  pat¬ 
tern  is  evolving  from  an  amorphous  triangular  region  at  1845 
MDT  to  a  bow  shaped  reflectivity  return.  The  strong  inflow 
winds  to  the  rear  of  the  bow  echoes  (south  for  Ml  and  east  for 
M2)  diminish  as  they  cross  the  reflectivity  gradient.  The  wind 
emerges  on  the  other  side  at  only  a  fraction  of  its  former  mag¬ 
nitude.  There  must  be  a  momentum  and  kinetic  energy  transfer 
downward  in  the  downdrafts  to  account  for  this  discrepancy.  The 
vertical  velocity  pattern  (Fig.  6.37),  shows  that  the  bow 
echoes  are  elongated  areas  of  downward  motion.  The  gust  front 
locations  at  lower  levels  remain  areas  of  strong  upward  motion. 
The  raicrobursts  are  surrounded  by  areas  of  upward  motion  caused 
by  the  influence  of  lower  level  gust  fronts  and  other  conver- 
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ger.ce  zones.  The  magnitude  of  the  updrafts  and  downdrafts  has 


increased  from  the  values  present  at  1845  MDT.  The  perturbation 


pressure  at  1.25  km  (Fig.  6.38)  reveals  that  the  misocyclones 


are  accompanied  by  low  pressure  and  M3  is  also  an  area  of  low 


pressure,  even  though  the  cyclonic  misocyclone  is  not  present. 


Again,  since  M3  is  just  forming,  the  wind  field  is  unable  to 


adjust  to  the  mass  field.  Over  a  span  of  several  minutes,  M3  is 


likely  to  develop  a  misocyclone  circulation  as  Ml  and  M2  have. 


The  dramatic  drop  in  relative  perturbation  pressure  and  exten¬ 


sive  coverage  of  the  low  pressure  region  suggests  that  M2  may 


be  intensifying  the  fastest,  which  is  also  reflected  in  the 


increased  organization  of  M2.  The  1.25  km  temperature  perturba¬ 


tion  pattern  (Fig.  6.39)  suggests  that  warm  environmental  air 


is  channeled  in  a  south  to  north  flow  toward  the  Ml  microburst. 


The  center  of  the  misocyclone  is  warm-cored  as  the  entrained 


inflow  penetrates  to  the  center  of  the  misocyclone.  The  core  is 


surrounded  by  cold  air.  This  suggests  that  the  precipitation 


cooled  air  mixes  with  the  warm  dry  environmental  air  to  produce 


an  enhanced  area  of  cold  air  ringing  the  misocyclone  at  its 


mature  state  (Ml).  The  newly  organized  misocyclone  (M2)  appears 


to  have  a  bipolar  temperature  field.  One  warm,  associated  with 


the  environmental  air,  and  one  cold,  associated  with  precipita¬ 


tion  cooling  effects,  are  noted.  No  evidence  of  misocyclone 


occlusion  by  cold  air  wrap  around  is  seen.  Ml  is  thought  to  be 
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close  to  total  occlusion  and  an  older  entity  than  Ml  or  M3. 
The  newly  formed  diffuse  microburst,  M3,  does  not  have  a  well- 
defined  misocyclone  at  1850  MDT,  but  one  appears  to  be  forming. 
The  diffuse  circulation  is  cold  cored,  and  is  associated  with  a 
descending  core  of  precipitation  and  precipitation  cooled  air. 
The  thermodynamic  quantities  for  the  tendency  case,  reveals 
that  difference  consists  of  a  greater  pressure  gradient  in  the 
tendency  case  (0.2  rab/km  vs.  0.25  mb/km)  for  a  comparable  spa¬ 
tial  distance  across  Mi  (Fig.  6.40).  The  perturbation  tempera¬ 
ture  pattern  (Fig.  6.41)  features  a  similar  pattern  as  featured 
in  the  non-tendency  case.  The  plan  view  may  be  complimented  by 
cross-sections  through  areas  of  interest  to  provide  a  three 
dimensional  view  of  the  storm  structure. 

6.2  Cross-sections  for  the  Complex  Cases 

The  complex  cases  require  further  kinematic  analysis  to 
determine  the  extent  of  the  interactions  between  the  micro¬ 
bursts.  Two  cross-sections,  one  from  northwest  to  southeast 
(NW-SE)  and  a  second  from  west  to  east  (W-E),  were  taken  across 
the  domain,  see  lines  AB  and  CD  in  Fig.  6.13.  The  following 
subsections  will  examine  these  cross-sections  at  times  1845  MDT 


and  18?  '  MDT,  5  August  1982. 
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Fig.  6.40  Same  as  Fig.  6.28,  except  for  1.25  km  with  ten¬ 
dency  included,  and  for  a  0.50  km  grid  spacing. 
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5  . £ .  !_  Cross-sect  ion  NV-SJ!  for  184-5  MPT ,  5.  August  1982 

The  1845  MDT  wind  flow  field  (Fig.  6.42)  shows  that  the 
fast  flow  is  generally  from  the  right  (southeast)  to  the  left 
(northwest)  with  an  area  between  (0,-26)  and  (-2,-24),  where 
the  flow  forms  a  cat  eye  type  turbulent  circulation.  It  is  a 
vortex  ring  resulted  by  the  interaction  between  the  microburst 
outflow,  and  the  environmental  inflow  (e.g.,  Fujita,  1988;  Lin 
and  Hughes,  1987). This  circulation  is  the  cross-section  of  the 
northern  outflow  boundary  of  Ml.  The  horizontally  spreading 
outflow  from  the  raicroburst  acts  as  an  obstacle  to  the  flow, 
but  the  flow  induced  a  secondary  circulation  on  the  microburst 
similar  to  what  occurs  in  separation  as  flow  moves  across  a 
valley  from  an  elevated  ridge.  Notice  that  at  (-6,-19)  some 
evidence  of  downward  motion  exists  at  1.25  km.  The  right  side 
of  Ml  is  where  the  gust  front  occurs.  The  convergence  is  ampli¬ 
fied  at  the  surface  by  the  return  flow  with  the  secondary  cir¬ 
culation.  To  the  left  of  Ml,  the  flow  descends  and  the  secon¬ 
dary  circulation  accentuates  the  downward  motion.  Perhaps  if 
this  circulation  extends  through  a  deep  enough  layer,  it  may 
induce  convection  on  the  right  side  through  focused  upward 
motion  and  further  raicroburst  development  on  the  left  side  by 
focused  downward  motion.  The  vertical  velocity  field  confirms 
upward  motion  (Fig.  6.43)  with  the  gust  front  on  the  right  side 
of  the  raicroburst  at  (-0.5,  -24.5),  and  downward  motion  on  the 
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ross-section 


left  side  at  (-2,-23).  Another  area  of  weak  downward  motion 
exists  aloft  at  (-6,-19),  which  is  associated  with  the  descend¬ 
ing  and  expanding  microburst  M2.  The  perturbation  pressure  gra¬ 
dient  (Fig.  6.44)  reveals  that  an  area  of  highest  pressure  on 
the  right  side  of  the  cross-section  accelerates  the  flow  to  the 
left,  through  the  strong  perturbation  pressure  gradient  between 
the  perturbation  pressure  high  and  the  low  associated  with  the 
gust  front.  Notice  that  once  the  flow  crosses  the  microburst 
outflow  it  descends  and  mass  accumulates  at  (-1 . 5 , -23 . 3 ) .  This 
creates  a  maximum  value  of  perturbation  pressure,  which  in  turn 
creates  a  reverse  pressure  gradient;  therefore,  the  flow  rev¬ 
erses  from  left  to  right  at  low  levels  toward  the  pressure 
deficit  associated  with  the  gust  front.  The  temperature  pattern 
(Fig.  6.45)  also  reveals  some  interesting  features.  The  air  to 
the  right  of  (0,-25)  is  warm  environmental  air.  As  the  air 
enters  the  gust  front  it  is  forced  to  rapidly  rise  and  adiabat- 
ically  cool  to  form  a  cold  pocket  up  to  0.6  km.  As  the  flow 
becomes  more  horizontal  west  of  the  most  intense  cooling,  the 
cooling  is  reduced  and  the  cold  pocket  is  less  intense.  Above 
0.5  km  at  (-2.5,  -22.5),  a  warm  pocket  in  conjunction  with  the 
descending  air.  As  the  air  nears  the  ground,  the  flow  becomes 
more  horizontal  and  the  downward  velocity  weakens  until  a  cool 
anomaly  is  formed  near  the  surface,  probably  due  to  an  evapora¬ 
tive  cooling  of  the  descending  air.  These  patterns  will  be 


6-2-2  Cross-section  NW-SS  for  1850  MPT,  5  August  1982 

The  1850  MDT  wind  flow  field  (Fig.  6.46)  shows  that  the 
cat  eye  pattern  associated  with  Ml  is  intact  with  a  slightly 
more  organized  circulation.  A  rotor  aloft  appears  to  be  forming 
at  (-4,-22)  at  1.0  km.  This  feature  forms  in  conjunction  with 
the  descent  of  M2  at  (-6,-20).  Note  that  the  microburst’s 
effect  on  the  wind  field  does  not  extend  to  the  lowest  level. 
Winds  at  0.25  km  are  still  dominated  by  a  right  to  left  com¬ 
ponent,  although  winds  to  the  southeast  of  M2  at  upper  levels 
do  show  convergence  at  (-5,-21)  and  a  left  to  right  flow.  The 
vertical  velocity  field  (see  Fig.  6.47)  shows  the  result  of  the 
intensification  and  organization  of  Ml  and  M2.  The  gust  front 
again  shows  a  +2  m/s  value  of  upward  motion.  On  the  other  side 
of  the  cat  eye  circulation,  a  value  of  -2  m/s  occurs.  At  (-4,- 
22)  the  convergence  zone  between  Ml  and  M2  produces  a  zone  of 
upward  motion,  mostly  aloft  and  above  0.5  km.  M2  shows  a 
dramatic  increase  in  vertical  velocity  to  -3  m/s  at  1.25  km. 
The  pressure  pattern  (Fig.  6.48)  shows  that  the  gust  front 
associated  with  Ml  is  shown  as  a  weak  trough.  The  high  pressure 
associated  with  Ml  only  extends  to  0.5  or  0.75  km.  M2  by  con- 
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Fig.  6.46  Same  as  Fig.  6.42,  except  for  1850  MDT. 
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;rast  extends  through  the  entire  layer,  reflecting  the  strong 


downward  motion  and  accumulation  of  mass  through  the  entire 


boundary  layer.  Between  Ml  and  M2  a  fairly  sharp  trough  exists. 


The  pressure  pattern  would  predict  a  convergence  between  Ml  and 


M2  at  all  levels,  but  the  convergence  does  not  extend  to  the 


lowest  level.  A  possible  explanation  could  be  the  dominance  of 


the  terrain  induced  flow,  which  is  primarily  felt  at  the  lowest 


level  for  this  terrain.  A  strong  gradient  of  pressure  to  the 


left  of  M2  accelerates  the  flow  to  that  direction.  The  tempera¬ 


ture  field  (see  Fig.  6.49)  also  shows  the  changes  in  structure 


of  Ml  and  M2.  The  cool  perturbation  temperature,  associated 


with  the  gust  front  near  Ml,  is  still  present.  The  convergence 


zone  above  0.5  km  and  between  Ml  and  M2,  reveals  a  cool  pocket, 


as  the  rising  aijr  cools  by  adiabatic  expansion.  A  warm  anomaly 


occurs  in  the  downward  flowing  air  in  the  Ml  cat  eye  circula¬ 


tion.  Evidently,  the  warmer  air  that  was  aloft  at  1845  MDT  has 


reached  the  surface  and  the  increased  downward  motion  enhances 


corapressional  warming.  M2  is  associated  with  an  approximate 


neutral  temperature  anomaly  aloft,  but  a  negative  anomaly 


occurs  at  lower  levels.  An  explanation  for  this  phenomenon  may 


be  that  as  the  air  descends  from  aloft,  it  mixes  at  the  lowest 


level,  where  environmental  flow  and  terrain  flow  dominate.  Mix¬ 


ing  at  the  lowest  level  ensues  as  the  moist  downburst  air  mixes 


with  the  environmental  air  at  the  0.25  to  0.50  km  level.  The 
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results  with  the  tendency  included  causes  a  noticeable  inten¬ 
sification  of  the  pressure  gradient  of  Ml  versus  M2  (Fig. 
6.50).  M2  is  approximately  the  same  in  magnitude  and  dimen¬ 
sions.  The  temperature  field  (Fig.  6.51)  reveals  that  the  cold 
anomaly  associated  with  the  upward  motion  with  the  convergence 
zone  between  Ml  and  M2  is  intensified  by  inclusion  of  tendency. 
The  warm  anomaly  associated  with  Ml  is  reduced  and  the  cold 
anomaly  associated  with  M2  is  also  reduced  in  magnitude  but  not 
in  sign.  An  east  to  west  cross-section  passing  through  Ml  and 
M3  will  be  examined  next  to  investigate  the  interaction  between 
those  raicrobursts. 


6 . 3_.  3.  Cross-section  W-E  for  1845  MPT ,  5  August  1982 

The  wind  flow  pattern  at  1845  MDT  (Fig.  6.52)  for  the  W-E 
cross-section  shows  a  right  (east)  to  left  (west)  flow  with  an 
area  of  strong  downward  motion  at  x  -  -2.  The  downward  motion 
is  associated  with  the  Ml  raicroburst.  Not-.-  that  strong  downward 
motion  extends  to  the  top  of  the  layer.  Also  note  that  the 
environment  1  air  mixes  with,  and  becomes  part  of,  the  down- 
draft  from  the  right  side  of  the  downdraft.  The  warm  dry 
environmental  air  contacts  the  precipitation  loaded  downdraft, 
and  undoubtably  precipitation  cooling  is  the  result.  The  down¬ 


ward  rushing  air  is  sufficiently  slow  to  allow  environmental 


Fig.  6.52  W-E  cross  -sect io 
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0.  Notice  that  the 


tively  warm  air  extends  from  x  =  4 
prevailing  right  to  left  flow  advects  the  warm  environmental 
air  into  the  finger  like  projection  at  x  -  0  at  a  height  of 
0.75  to  1.0  kra.  A  cold  anomaly  occurs  with  the  microburst  down- 
draft,  which  is  especially  noticeable  near  the  0.25  km  level.  A 
much  weaker  positive  anomaly  exists  in  the  warmer  environmental 
air  west  of  the  Ml  microburst.  This  feature  apparently  is 
caused  by  compressionally  warmed,  downwardly  induced  air  flow 
to  the  lee  of  the  microburst.  Obstacle  flow  is  probably  the 
cause.  Another  cold  anomaly  exists  at  the  west  end  of  the  grid. 
The  origin  of  this  feature  is  not  known.  It  may  be  the  remains 
of  a  previous  microburst  event.  The  second  time  period  will  be 
used  to  document  changes  in  the  vertical  structure  of  the 
microbursts. 

6.3.U  Cross-section  W-E  for  1850  MPT,  5^  August  1982 

The  wind  flow  at  1850  MDT  (Fig.  6.56)  is  characterized  by 
an  increase  in  turbulent  motion,  as  the  microbursts  intensify 


and  become  better  defined.  The  main  microburst,  Ml,  appears  to 
be  split  into  two  areas  of  strong  downward  motion  on  either 


side  of  x  -  0.  This  twinning  may  be  due  to  an  accumulation  of 
mass  at  lower  levels  having  a  higher  density,  because  mixing 
with  warm  dry  environmental  air  will  cool  as  precipitation  is 


.56  Same  as  Fig.  6.52,  except  for  1850  MDT. 


mixed  with  it.  The  denser  air  then  accumulates  at  the  lowest 
levels.  As  warmer  lighter  air  descends  from  aloft,  it  flows 
down  the  flanks  of  the  pocket  of  denser  air,  leaving  a  rela¬ 
tively  undisturbed  region  at  the  microburst  center.  The  verti¬ 
cal  velocity  field,  shown  in  Fig.  6.57,  also  illustrates  some 
major  changes.  At  x  *  4  an  incipient  downdraft  is  observed  at 
levels  above  0.50  km.  This  feature  is  part  of  the  M3  micro¬ 
burst.  Another  area  of  weak  downward  motion  occurs  above  0.75 
km  at  x  -  -6.  This  area  corresponds  with  the  outermost  fringe 
of  the  M2  microburst  system.  It  is  likely  associated  with  the 
downward  and  inwardly  flowing  air  entrained  by  the  misocyclone 
associated  with  M2.  At  1845  MDT  the  misocyclone  was  very  disor¬ 
ganized  and  this  feature  was  not  present.  At  1850  MDT  it  became 
better  organized  and  larger  in  area.  The  vertical  velocity  dep¬ 
icts,  from  right  to  left;  (1)  an  area  of  weak  descending  air 
associated  with  M4,  (2)  a  narrow  corridor  of  rising  motion 
associated  with  the  gust  front,  (3)  a  double  structure  micro- 
burst  downdraft,  which  has  increased  in  magnitude  from  2  to  3 
m/s  .showing  Ml  and  M3,  (4)  an  area  of  weak  upward  motion  asso¬ 
ciated  with  the  gyre  between  Ml  and  M2  and  (5)  the  weak  area  of 
downward  motion  between  x  -  -4  and  -6  associated  with  the 
fringe  of  M2.  The  descent  of  M2  and  intensification  of  Ml  has 
produced  a  gyre  between  them  aloft  at  x  -  -3.  A  similar  event 
occurred  during  the  evolution  of  the  Ml  cat  eye  circulation  and 
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M2  nicroburst  for  the  NV-SE  cross-section.  Apparently,  the 
gyres  and  associated  strong  turbulence  develop  and  maintain 
themselves  when  directional  shear  is  established  between  two 
microburst  circulations.  At  least  for  weak  microbursts,  an 


interaction  between  raicrobursts  may  be  useful  in  creating  gyres 
and  strong  turbulence.  The  strongest  turbulence  can  thus  be 
expected  to  occur  between  microburst  events  at  the  0.75  and 
1.25  km  level.  The  pressure  field  (Fig.  6.58)  shows  a  trough  at 
x  -  2  extending  upwards.  The  centers  of  Ml  and  M2  are  located 
at  x  -  -1.5  and  x  -  1.5,  respectively.  The  microburst  extends 


upwards  to  0.75  km  as  a  distinct  system.  The  temperature  field 
(Fig.  6.59)  shows  a  sharp  thermal  trough  at  x  »  3 ,  associated 
with  warm  air  in  the  updraft  region.  Above  0.75  km  at  x  -  1, 


the  temperature  anomaly  is  warm.  This  suggests  that  the  des¬ 
cending  downdraft  is  initially  warm,  but  at  approximately  1  km 
dry,  warm  environmental  air  is  able  to  enter  and  mix  with  the 


descending  unsaturated  air.  The  evaporatively  cooled  air 
becomes  increasingly  cool  as  it  descends.  The  pressure  pattern 
with  tendency  (Fig.  6.60)  suggests  that  M2  is  amplified  and 
better  defined.  The  gyre  between  Ml  and  M2  becomes  defined  as  a 
trough  of  low  pressure  and  the  amplitude  of  the  Ml  gust  front 
to  the  east  of  Ml  is  less  sharply  defined.  The  temperature  pat¬ 
tern  (Fig.  6.61)  confirms  the  features  present  in  the  non¬ 
tendency  case,  and  it  also  increases  the  magnitude  and  coverage 
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of  the  told  anomaly  associated  with  M2. 

The  three-dimensional  view  of  the  microburst  for  two  dif¬ 
ferent  time  periods  reveals  an  intensifying  microburst  struc¬ 
ture  consisting  of  a  descending  reflectivity  core,  accompanied 
by  increasing  downward  motion.  The  downward  motion  is  fed  by 
the  entrainment  of  dry,  warm  environmental  air  from  the  upper 
layers.  This  high  speed  air  is  fed  into  the  bow  echo  and  much 
of  the  momentum  and  kinetic  energy  appears  to  be  transferred 
downward  in  the  downdrafts.  The  initially  warm  downdraft  cores 
are  ventilated  by  the  inwardly  spiraling  environmental  air  fed 
in  by  the  misocyclone.  The  misocyclone  intensifies  and  becomes 
vertically  compressed  as  the  vertical  velocity  increases,  as 
the  vorticity  of  the  column  becomes  stretched.  The  stretched 
column  spins  down  when  the  downdraft  nears  the  surface.  The 
cyclonic  circulation  reaches  a  neutral  point  near  the  0.75  km 
level  then  develops  an  anticyclonic  circulation  beneath  that 
point.  The  5  August  1982  case  contains  hail  and  this  may  be 
critical  to  initiating  the  raicroburst  in  the  relatively  moist 
environment  of  this  storm.  The  14  July  1982  case  may  not  have 
any  significant  amount  of  hail.  This  may  be  expected,  since  it 
forms  in  a  drier  environment;  hence,  evaporative  cooling  alone 
can  initiate  and  maintain  the  microburst.  The  simple  case  does 
not  have  a  well-defined  misocyclone.  This  may  be  due  to  the 
fact  that  the  stronger  horizontal  winds  results  in  a 


modification  of  the  misocyclone  structure.  Both  the  simple  and 
the  complex  cases  feature  entrainment  of  dry,  warm  high-speed 
air  into  the  windward  side  of  the  bow  echo.  The  interactions 
between  the  microburst  and  its  environment  are  limited  to  gust 
fronts  in  the  simple  case.  The  complex  case  interactions 
results  in  microbursts  that  enhance  gust  front  developments, 
development  of  enhanced  areas  of  upward  motion  rotors  and  a 
microburst  embedded  in  a  stronger  microburst.  These  observa¬ 
tions  may  be  further  justified  by  examining  the  results  of  the 
eddy  kinetic  energy  budgets  and  momentum  budgets. 

6.4  The  Momentum  Budgets 

The  momentum  budgets  presented  in  this  section  are  for  the 
1845  and  1850  MDT,  5  August  1982  cases.  In  addition  some  com¬ 
parisons  will  be  made  with  the  1649  MDT  14  July  1982  case,  but 
the  results  will  not  be  shown  because  of  the  relative  simpli¬ 
city  of  the  result  and  because  of  the  emphasis  placed  on  the 
complex  case  for  this  research.  The  cases  will  be  examined  in 
two  parts.  The  first  part  will  consist  of  the  local  turbulent 
eddy  momentum  fluxes.  The  second  part  will  consist  of  the  total 
momentum  flux  budget,  as  patterned  after  LeMone  (1983). 


6.4. _1  The  Local  Turbulent  Eddy  Momentum  Fluxes 

The  vertical  components  of  local  turbulent  u-  and  v- 
momentum  fluxes  (u'w'  and  v'w')  produce  a  picture  of  how  momen¬ 
tum  is  being  transferred  by  the  turbulent  eddies.  The  following 
analysis  will  show  how  horizontal  eddy  momentum  is  being 
transferred  by  vertical  motion,  which  features  are  predom¬ 
inately  responsible  for  the  momentum  transfer,  and  how  features 
are  temporally  evolving. 


The  Local  Turbulent  Momentum  Flux  at  1845  MPT 

The  first  eddy  flux  quantity  examined  will  be  the  u'w' 
momentum  flux,  along  the  vertical  cross-section  (see  Fig. 
6.13).  A  gain,  represented  by  a  positive  anomaly  centered  at 
0.75  km,  near  Ml  (Fig.  6.62),  is  caused  by  downward  motion  (w' 
<  0)  of  westward  (u'  <  0)  momentum  air.  General  gains  occur  in 
a  column  of  air  above  microburst  Ml  for  the  same  reason.  To  the 
right  and  above  the  gust  front  (dashed  line),  a  small  intense 
area  of  loss  of  momentum  occurs.  This  area  corresponds  to  west¬ 
ward  moving  air  within  the  updraft,  showing  the  u-component  of 
eddy  momentum  is  being  transferred  downward  (u'w’  <  0)  in  the 
updraft  column.  The  momentum  transfer  appears  to  be  one  of  loss 
(downward)  from  the  updraft  region  with  gain  (upward)  in  the  Ml 


microburst  region.  Recall  from  the  kinematic  analysis  of  the 


west  to  east  cross-sections,  the  area  between  the  loss  and  gain 
centers  is  the  area  of  westward  horizontal  flow  from  the  loss 
region  to  the  gain  region.  The  areas  of  positive  u'w'  show  a 
discrete  maximum  within  the  atmospheric  column  over  Ml.  The  M2 
microburst  is  not  very  active  in  terms  of  momentum  flux  at  1845 
MDT.  There  is  a  weak  loss  of  momentum  near  the  top  of  the 
domain  at  1.25  km,  due  to  the  downward  motion  of  eastward  mov¬ 
ing  air. 

The  v'w'  momentum  flux  (Fig.  6.63)  reveals  that  the  gust 
front  is  the  major  sink  (downward  transport)  of  v-momentum,  as 
it  was  for  u-moraentum.  Notice  that  the  sink  is  at  0.5  km 
instead  of  0.75  km  for  the  u'w'  case.  Since  the  gust  front  is 
mostly  a  discontinuity  in  the  v-component  of  the  wind  for  this 
case,  and  the  discontinuity  only  extends  to  between  0.50  and 
0.75  km,  the  maximum  sink  occurs  at  0.5  km.  Also  note  that  a 
negative  value  of  u'w'  above  0.75  km  corresponds  with  a  weak 
positive  value  of  v'w'.  This  implies  some  u-momentum  may  be 
converted  to  v-moraentum  in  this  area.  Notice  that  within  the  Ml 
atmospheric  column,  a  positive  value  of  v'w'  occurs  at  the 
lowest  levels  with  negative  values  above.  This  suggests  that 
v-raomentura  is  lost  in  the  vicinity  of  the  misocyclone,  and  is 
transported  downward  to  the  lower  boundary  layer,  where  the  Ml 
microburst  is  located.  The  M2  circulation  shows  a  source 
(upward  transport)  of  v'w'  momentum  flux  at  the  highest  level, 
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corresponding  to  the  southward  v-momentum  (v1  <  0)  being 
transferred  downward  within  the  downdraft  column  of  M2.  The 
kinematic  analysis  showed  a  general  intensification  of  the 
microburst  features  by  1850  MDT.  An  analysis  of  the  turbulent 
momentum  fluxes  should  provide  information  on  how  that 
transformation  occurred. 


6 . 4. 1_.  2^  The  Local  Turbulent  Momentum  Flux  at  1850  MDT 

The  u'w*  momentum  flux  for  1850  MDT  (Fig.  6.64)  shows  that 
the  gust  front  remains  active,  but  it  is  no  longer  the  area  of 
greatest  momentum  flux.  The  atmospheric  column  above  M2 
represents  a  large  sink  (downward  transport)  of  u-momentum, 
indicative  of  the  eastward  u-momentum  (u'  >  0)  being 
transferred  downward  within  M2.  Most  of  the  loss  occurs  above 
0.75  km  with  the  largest  value  at  1.0  km.  The  explosive  growth 
of  the  sink  in  u'w'  occurs  simultaneously  with  the  rapid  organ¬ 
ization  and  growth  of  the  M2  misocyclone  at  those  levels.  The 
implication  of  this  development  suggests  that  the  downward 
vertical  motion  has  intensified  the  removal  of  mass  and  momen¬ 
tum  from  levels  above  0.75  km,  and  transported  it  towards  the 
surface.  Ml  has  likewise  intensified  with  a  large  positive  u- 
momentum  accumulation  near  0.75  km  with  a  weak  negative  sink  in 


u-momentum  near  1.25  km.  The  v'w'  (Fig.  6.65)  momentum  flux 


Fig.  6.64  Same  as  Fig.  6.62,  except  for  1850  MDT. 

Fig.  6.65  Same  as  Fig.  6.62,  except  for  1850  MDT  and  v'w 
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reveals  that  the  large  sink  of  u'w'  associated  with  M2  is  coin¬ 
cident  with  the  large  source  of  v'w'.  This  implies  a  conversion 
of  momentum  from  u'w'  to  v'w'  by  the  misocyclone  microburst 
circulation.  The  momentum  flux  v'w'  for  the  Ml  microburst  seems 
static  in  comparison  to  the  explosive  growth  exhibited  by  M2. 
It  appears  that  Ml  is  more  nearly  quasi-stationary  than  M2,  and 
likely  near  the  peak  of  its  development  process. 

A  pattern  of  general  upward  or  downward  flux  of  momentum 
for  the  entire  domain  may  be  seen  by  examining  the  vertical 
profiles  for  1845  and  1850  MDT  (Fig.  6.66).  The  <u'w'>  profile 
at  1845  MDT  suggests  a  flux  of  momentum  upward  to  misocyclone 
levels  and  downward  to  the  microburst  level.  This  result  sug¬ 
gests  that  the  misocyclone  is  losing  momentum,  the  microburst 
is  gaining  momentum,  and  the  mid-level  is  largely  losing  momen¬ 
tum  with  little  gain  to  offset  the  loss.  This  result  shows  that 
the  misocyclone  circulation  provides  momentum  to  the  low  lev¬ 
els.  The  misocyclone  circulation  is  also  the  site  of  accumula¬ 
tion,  as  momentum  is  transferred  from  below  in  the  updrafts 
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which  is  horizontally  spiraled  inward  into  the  misocyclone  cir¬ 
culation.  At  1850  MDT  the  pattern  of  downward  flux  to  low  lev¬ 
els  continues  with  a  slight  reversal  at  0.25  km.  The  v'w'  pro¬ 
file  at  1845  MDT  suggests  a  pattern  of  momentum  flux  to  higher 
and  lower  levels,  but  a  general  transfer  downward  is  evident. 
The  pattern  at  1850  MDT  changes  dramatically  as  Ml  intensifies. 
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The  slope  shows  large  losses  aloft  with  the  smallest  or  loss  at 
0.25  km.  The  pattern  coincide  with  the  kinematic  analysis.  The 
kinematic  analysis  showed  that  the  0.75  km  level  was  virtually 
devoid  of  eddy  motion,  which  is  predominately  composed  of  the 
v-coraponent  of  motion;  therefore,  u-momentum  is  predominately 
present  in  the  turbulent  flow  above  and  below  0.75  km,  while 
u-momentum  is  present  at  all  levels  with  maximum  values  present 
in  the  strong  winds  at  0.25  and  1.25  km.  U-momentum  is  continu¬ 
ously  extracted  from  levels  above  0.50  km  and  is  added  to  the 
misocyclone  for  transport  downward  to  enhance  the  eddy  motion 
of  the  microburst. 

While  the  local  turbulent  momentum  budgets  are  useful  in 
tracking  the  sinks  and  sources  of  the  eddy  momentum,  a  momentum 
budget  containing  representations  of  the  underlying  physical 
mechanisms  is  necessary  to  determine  the  source  of  the  momentum 
transformations . 


6.4.2  The  Total  Momentum  Budgets 

The  momentum  budgets,  as  derived  by  LeMone  (1983),  will  be 
used  to  examine  the  large  scale  domain  given  by  the  13  by  13 
array  within  the  15  X  15  km  array.  A  second  smaller  domain, 
taken  to  isolate  the  momentum  budget  processes  within  Ml  was 


taken  across  a  4  X  4  km  array  centered  on  the  Ml  microburst  for 


the  5  August  1982  case,  and  the  singular  microburst  for  the  14 
July  1982  case.  This  nested  approach  may  help  isolate  mechan¬ 
isms  associated  with  the  interactions  of  the  microbursts  versus 
those  mechanisms  that  are  important  to  only  the  microbursts. 
The  1845  and  1850  MDT,  5  August  1982  cases  will  be  examined 
sequentially  with  comparisons  with  the  1649  MDT,  14  July  1982 
case,  which  will  not  be  shown. 

6 .4 . 2  ._1  The  Momentum  Budgets  for  the  1845  MDT  Case 

The  large  domain,  momentum  budgets  at  1845  MDT  (Tables 
2,3)  reveal  that  the  balance  in  the  momentum  budget  is  largely 
maintained  by  a  balance  between  the  buoyancy  term  (B),  which 
acts  as  the  source  and  the  pressure  velocity  interaction  term 
(P)  that  acts  as  a  sink.  Since  a  large  amount  of  the  outflow 
from  the  cold-cored  microbursts,  particularly  Ml,  flows  to  the 
west,  the  buoyancy  term  is  positive  at  low  levels.  At  misocy- 
clone  levels  much  of  the  cold  inflow  is  associated  with  east¬ 
ward  moving  air;  therefore,  the  buoyancy  becomes  a  negative, 
for  the  v-raomentura  budget. A  different  scenario  occurs  in  the 
v-momentum  case  with  the  large  warm  anomaly  north  of  Ml  associ¬ 
ated  with  northward  moving  air  at  low  levels.  The  pressure 
velocity  interaction  term  is  a  sink  for  both  the  u  and  v- 
momentum  budgets  at  the  lowest  two  levels.  The  u-momentura 
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budget  accurately  reflects,  that  at  the  outer  edges  of  the 
microburst,  upward  motion  is  expected.  This  does  not  hold  above 
0.50  km,  because  the  microburst  no  longer  dominates.  The 
environmental  wind,  which  is  predominately  northerly  in  direc¬ 
tion,  is  sufficiently  strong  to  modify  the  misocyclone  circula¬ 
tion  in  favor  of  a  positive  contribution  to  the  momentum  budget 
from  the  pressure  wind  interaction  term. 

The  shear  (S)  and  transport  (Tr)  terras  are  secondary  con¬ 
tributors  to  the  u-and  v-momentum  budgets  at  low  levels.  The 
transport  term  features  a  source  in  the  layer  between  the  sur¬ 
face  and  0.5  km,  and  a  sink  aloft  for  both  cases.  The  shear 
term  reveals  a  source  of  momentum  at  all  levels  except  at  the 
highest  level.  The  misocyclone  circulation  is  extracting  mean 
u-momentum  from  aloft  and  transporting  it  downward.  The  v- 
momentum  shear  term  is  increasingly  a  source  of  momentum  with 
height.  This  suggests  that  the  shear  term  transports  u- 
raoraentum,  extracted  from  the  mean  flow,  downward  in  the 
microburst/raisocyclone  circulation.  Since  the  environmental 
flow  is  south  to  north  and  the  environmental  component 
increases  with  height,  the  vertical  motions  associated  with  the 
microburst  are  able  to  add  momentum  to  the  eddy  scale  from  the 
mean  scale  to  help  support  the  v-raomentum  budget. 

The  tertiary  contributors  to  the  momentum  budgets  are  the 
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horizontal  divergence  of  u-  and  v-transport  and  the  vertical 
divergence  of  w-transport.  The  most  significant  aspect  is  that 
the  divergence  of  u-momentura  contributes  more  positively  to  v- 
momentum  than  the  u-momentum  budget,  and  is  a  sink  for  the  v- 
momentum  budget.  This  suggests  there  is  an  exchange  or  conver¬ 
sion  of  momentum  between  the  u  and  v  direction.  The  larger 
scale  interactions  reveal  that  buoyancy  balances  the  pressure 
velocity  interactions  as  the  main  contributors.  The  secondary 
contributors  of  shear  and  transport  reveals  that  the  misocy- 
clone  extracts  anargy  from  the  mean  u-wind  flow  through  the 
shear  term,  with  a  gain  in  v-momentum  occurring  at  all  levels 
as  a  result  of  the  vertical  motion  of  the  updrafts  and  down- 
drafts  acting  on  a  negative  mean  v-wind  shear.  The  transport 
terra  confirms  the  results  from  the  local  momentum  flux.  The 
momentum,  both  u  and  v,  is  being  extracted  by  the  misocyclone 
at  upper  levels  and  transported  downward  to  the  raicroburst. 
What  is  interesting  to  note,  is  that  the  u-momentum  experiences 
a  loss  of  momentum  at  all  levels,  while  the  v-momentum  budget 
experiences  a  gain  at  all  levels  but  the  lowest.  The  micro¬ 
burst  influences  the  domain  of  interest,  but  to  examine  only 
the  microburst  structure  a  smaller  domain  is  required.  By  exa¬ 
mining  the  microburst  structure  it  may  be  possible  to  pinpoint 
the  relative  importance  of  transport  in  the  overall  microburst 
structure . 


The  Ml  microburst  momentum  budgets  (Tables  4,5)  reveals 
that  the  buoyancy  and  velocity  pressure  interaction  terms 
remain  important,  but  the  shear  and  transport  terms  become  pri¬ 
mary  contributors  at  the  microburst  levels.  Even  the  Dx  and  Dy 
terms  become  significant  at  the  misocyclone  level.  Another 
development  of  note  is  the  emergence  of  the  local  time  deriva¬ 
tive  term,  as  a  primary  term.  This  confirms  that  the  most  rapid 
changes  in  the  larger  scale  domain  are  a  result  of  changes  in 
the  microbursts.  The  microbursts  largely  dominate  the  larger 
scale  flow  for  this  meso- /scale  domain. 

The  buoyancy  term  remains  a  source  for  the  v-momentum 
case.  It  becomes  very  large  below  0.50  km,  because  the  cold 
temperature  anomaly  is  combined  with  the  strongest  southward 
motion.  The  buoyancy  term  for  the  u-momentum  case  becomes  a 
negative  at  upper  levels,  as  the  warm  anomalies  become  combined 
with  easterly  flow,  and  the  entrainment  flow  within  the  misocy¬ 
clone  is  predominately  easterly.  At  larger  scales  the  broader 
westerly  inflow  dominates.  At  the  smaller  scale  the  constricted 
inflow  near  the  misocyclone  center  is  easterly  and  a  negative 
contribution  is  made.  The  pressure  velocity  interaction  term 
behaves  in  a  manner  similar  to  the  large  scale.  The  shear  term 
becomes  a  more  important  component.  It  becomes  the  most  impor¬ 
tant  term  at  the  misocyclone  level  for  both  budgets.  Its  basic 
behavior  remains  unchanged,  but  its  importance  to  the  budget 
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increases,  particularly  at  the  misocyclone  level  where  the 


vertical  velocity  increases.  The  transport  term  also  becomes  a 
major  contributor,  particularly  at  the  misocyclone  level.  This 
development  confirms  the  role  of  the  misocyclone,  as  a  tran¬ 
sporter  of  momentum  downward  within  the  larger  circulation. 
This  implies  that  the  misocyclone  is  a  gatherer  and  transporter 
of  momentum  downward.  The  importance  of  Dx  and  Dy  in  the 
respective  budgets  of  the  u  and  v-momentum  budgets  at  the  miso¬ 
cyclone  level  suggests  that  a  large  amount  of  turbulent 
exchange  occurs  there.  More  importantly,  the  role  of  Dx  and  Dy 
in  the  respective  u  and  v-momentum  budgets  is  reversed  from  the 
large  scale  case.  The  Dx  term  becomes  a  source  for  the  u- 
momentum  budget,  and  Dy  becomes  a  source  for  the  v-momentum 
budget.  Apparently,  there  is  a  scale  interaction  between  the 
microburst/misocyclone  system,  and  the  larger  scale  environment 
within  this  terra.  The  environment  dominates  at  the  large  scale, 
where  a  momentum  exchange  between  the  u  and  v-momentum  budgets 
occur.  At  the  raicroburst  scale,  the  exchange  between  budgets 
does  not  appear  to  be  present  at  the  misocyclone  level,  but 
rather  a  direct  gain  or  loss  of  momentum  occurs  within  the  same 
directional  scale. 

The  results  of  this  examination  may  be  further  substan¬ 
tiated  by  examining  the  second  time  period,  1850  MDT. 


—  •  —  •'L-Il  The  Momentum  Budgets  for  the  1850  MPT  Case 

The  large  domain  budgets  at  1850  MDT  (Tables  6,7)  confirms 
the  results  from  the  1845  budgets.  The  buoyancy  term  balances 
the  pressure  velocity  interaction  term.  The  shear  and  transport 
terms  are  secondary  contributors  and  the  divergence  terms  (Dx, 
Dy  and  Dz)  are  tertiary  terms.  Even  small  scale  features,  such 
as  the  extraction  of  u-momentum  by  the  shear  term  at  microburst 
level,  remains  virtually  unchanged.  The  main  points  of  interest 
are  therefore  the  differences.  A  major  difference  is  the  weak¬ 
ening  of  the  magnitude,  but  not  the  sign  of  the  transport  term. 
This  suggests  that  either  the  intensification  of  the  micro¬ 
burst,  accomplished  by  enhanced  downward  transport  from  raisocy- 
clone  levels  may  be  leveling  off  or  a  counter  circulation  has 
developed.  The  main  intensification  phase  of  the  microbursts, 
particularly  M2,  may  be  over  by  1850  MDT.  M3  is  developing  at 
1850  MDT,  but  it  is  such  a  small  feature  that  Ml  and  M2  dom¬ 
inate  the  budget.  Another  difference  is  the  reversal  of  the 
role  that  Dx  and  Dy  play  in  the  momentum  budgets.  The  Dx  term 
now  contributes  positively  to  the  u-momentum  budget  and  Dy  con¬ 
tributes  positively  to  the  v-raoraentum  budget.  This  may  be  due 
to  the  expanding  and  intensifying  microbursts  dominating  these 
processes  at  the  expense  of  the  environment,  or  perhaps  a 
reversal  of  roles  played  by  the  microburst  and  the  larger  scale 
environment.  Only  an  examination  of  the  microburst  domain 
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budgets  may  provide  the  answer. 
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The  raicroburst  domain  budgets  at  1850  MDT  (Tables  8,9) 
reveal  that  the  Dx  and  Dy  term  retain  the  same  role  played  at 
1845  MDT.  Therefore,  the  growth  and  development  of  the  micro¬ 
burst  circulations  may  have  been  sufficient  to  dominate  the 
larger  scale  flow.  The  transport  term  increased  in  importance 
at  the  microburst  scale,  although  it  is  less  important  to  the 
largo  scale.  The  intensification  of  the  microbursts  have 
increased  the  transport  downward  within  the  microburst  domain; 
therefore,  the  likely  cause  of  the  weakening  in  the  magnitude 
of  the  transport  of  momentum  at  large  scales  is  due  to  feedback 
effects  brought  on  by  secondary  circulations.  As  the  micro- 
bursts  intensify,  the  gust  fronts  and  vertical  motions  between 
them  intensify.  This  was  discussed  in  the  kinematic  study  as 
the  intensification  and  creation  of  intense  gust  fronts,  verti¬ 
cal  rotors,  and  convergence  lines  brought  on  by  the  interac¬ 
tions  between  the  raicrobursts.  These  secondary  contributors 
counter  the  raicroburst  at  larger  scales,  by  transporting  some 
of  the  momentum  brought  down  from  aloft.  I  speculate  that  decay 
of  the  microburst/misocyclone  circulation  begins,  when  the 
counter  processes  cross  the  equilibrium  point  towards  a  posi¬ 
tive  flux  of  momentum  upward. 


In  comparison  to  the  5  August  1982  case,  the  1649  MDT,  14 


m 


V 
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July  1982  case  contains  the  following  characteristics.  Unlike 
the  5  August  1982  case,  the  balance  of  the  momentum  budget  is 
more  complicated.  The  balance  at  low  levels  is  maintained  by 
the  transport,  buoyancy,  and  pressure  velocity  tarns.  At  mid 
levels  the  balance  is  maintained  by  the  shear  and  buoyancy 
terms,  while  at  upper  levels  the  shear,  transport  and  buoyancy 
terms  dominate.  These  results  are  indicative  of  the  differences 
in  structure  between  the  5  August  and  the  14  July  1982  cases. 
The  latter  case  has  a  low  level  jet  at  3.5  km  (Lin  and  Hughes, 
1987),  a  downdraft  502  or  greater  in  magnitude,  and  an  updraft 
an  order  of  magnitude  larger  than  the  5  August  1982  case.  The 
result  is  that  the  larger  vertical  motion  and  strong  wind  shear 
provided  more  momentum  to  drive  the  conversion  of  momentum  from 
the  mean  to  the  eddy  flow  within  the  shear  term.  The  shear  term 
is  a  negative  for  the  u-raomentura  case,  and  a  positive  for  the 
v-momentura  case.  The  eddies  extract  energy  from  the  mean  v- 
momentura  and  add  it  to  the  mean  u-raoraentum.  The  transport  is  a 
maximum,  because  of  the  large  value  of  downward  vertical  velo¬ 
city  and  divergence  in  the  microburst  domain.  A  large  positive 
magnitude  at  the  low  levels  suggests  that  the  momentum  is  accu¬ 
mulating  below  0.50  km,  while  the  large  negative  magnitude  at 
the  upper  levels  reveals  that  momentum  is  transported  downward. 
The  pressure  velocity  term  is  a  minimum,  because  the  updraft 
and  downdraft  in  this  case  were  of  roughly  the  same  magnitude, 


and  they  tend  to  offset  each  other.  This  may  suggest  that  this 
microburst  is  entering  the  later  stages  of  its  development,  as 
the  secondary  circulations  assume  greater  importance.  The  com¬ 
parison  has  shown  that  microbursts  contain  unique  dynamics  that 
change  from  case  to  case.  A  study  of  many  cases  may  be  required 
to  definitively  identify  and  describe  the  complete  structure 
and  development  of  microbursts. 

The  research  shows  that  the  momentum  budgets  appear  to  be 
physically  consistent  with  the  kinematic  analysis.  There  is 
some  evidence  to  suggest  that  the  mean  flow  is  providing  sup¬ 
port  to  the  eddy  flow.  By  examining  the  eddy  kinetic  energy 
budget,  the  degree  of  the  support  and  its  physical  mechanisms 
are  the  objects  of  interest. 

6 . _5  The  Eddy  Kinetic  Energy  Budgets 

The  eddy  kinetic  energy  budget  will  be  examined  to  deter¬ 
mine  if  the  energy  transformations  agree  with  changes  in  the 
kinematic  and  momentum  budget  fields.  The  role  of  mean  and 
eddy  interactions  can  be  studied,  which  will  shed  light  on  how 
scale  interaction  influences  the  microbursts.  The  energy  budg¬ 
ets  first  examined  will  consist  of  the  vertical  flux  of  eddy 
kinetic  energy,  w'E'  .  This  terra  is  thought  to  be  important  in 
maintaining  the  microburst  through  transport  of  kinetic  energy. 
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E' ,  associated  with  the  gust  front,  remains.  The  structure  of 


the  Ml  microburst  is  perhaps  slightly  accentuated  by  an 
increase  in  the  area  coverage  of  the  positive  eddy  kinetic 
energy.  M2  by  contrast  reveals  a  lowering  and  intensification 
of  the  positive  anomaly  of  eddy  kinetic  energy.  The  vertical 
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flux  of  eddy  kinetic  energy  shows  that  the  growth  of  the  posi¬ 
tive  anomaly  in  E'  corresponds  to  the  explosive  growth  of  nega¬ 
tive  flux  of  E'.  The  kinematic  analysis  and  momentum  flux  budg¬ 
ets  examined  previously  suggested  that  M2  was  intensifying  more 


rapidly  than  Ml.  The  vertical  flux  of  E'  has  isolated  the 
mechanism  responsible  for  the  rapid  intensification  of  M2,  the 


raisocyclone.  The  rapid  intensification  of  M2  and  its  lowering 
by  vertical  contraction,  brought  on  by  spin  up  causes  a  rapid 
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increase  of  E'  and  a  lowering  of  the  feature.  The  intensified 
and  organized  misocyclone  is  able  to  efficiently  draw  in  dry 
environmental  air,  which  mixes  with  the  high  E'  air  within  the 
misocyclone.  The  added  negative  buoyancy  dramatically  increases 
the  downward  motion  with  a  dramatic  increase  in  W'E'.  The 
increased  downward  motion  stretches  the  column  and  spins  the 
misocyclone  at  a  faster  rate.  Undoubtably,  there  is  a  strong 
feedback  mechanism  between  the  misocyclone ' s  ability  to  entrain 
dry  air  and  increase  the  downward  motion,  and  the  role  that  a 
stronger  updraft  has  on  spinning  up  and  organizing  the  misocy¬ 
clone.  The  M2  misocyclone  was  weak  at  1845  MDT,  so  it  is  doubt- 
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ful  that  the  misocyclone  played  a  role  in  the  initiation  of  the 


microburst,  but  it  clearly  plays  a  role  in  the  intensification 


of  M2.  Ml  does  not  intensify  as  rapidly  as  M2  does.  Recall  that 


the  Ml  misocyclone  was  well-developed  at  1845  and  1850  MDT.  The 


WE'  profile  at  1850  MDT  looks  similar  to  the  1845  MDT  profile. 


The  1850  MDT  case  reveals  that  the  narrow  shaft  of  downward 


flux  of  E'  has  intensified  slightly.  Other  features  include  the 


increase  of  positive  area  in  WE'  between  Ml  and  M2.  As  M2 


increases  in  intensity,  it  becomes  better  able  to  oppose  the 


combined  outflow  of  Ml  and  M2,  since  M2  becomes  a  greater  obs¬ 


tacle  in  Ml's  outflow.  Another  interesting  feature  is  the 


decrease  in  W'E'  associated  with  the  gust  front.  This  may  be  a 


reflection  of  the  weakening  of  that  section  of  the  gust  front. 


The  E'  and  W'E'  profiles  were  useful  in  confirming  previous 


speculation  about  intensification  trends  in  the  Ml  and  M2 


raicrobursts.  It  also  confirmed  speculation  about  the  role  the 


misocyclone  circulation  plays  in  the  intensification  and 


maintenance  of  the  microburst  circulation.  It  does  not  identify 


specific  physical  processes  that  contribute  to  the  growth  of 


the  microbursts.  This  missing  key  may  confirm  the  role  of  dry 


entrainment  in  the  misocyclone  circulation.  The  missing  key  may 


be  provided  by  examining  the  eddy  kinetic  energy  budget  for  the 


1845  and  1850  MDT  time  periods. 
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6.^.2  The  Complete  Eddy  Kinet ic  Energy  Budget 

The  eddy  kinetic  energy  budget  will  be  examined  at  the 
1845  and  1850  MDT  time  periods  and  for  the  1649  MDT  case  too. 
At  1850  MDT  two  cases  will  be  considered.  One  case  includes 
friction,  and  one  case  does  not.  The  eddy  kinetic  energy  budget 
will  be  further  subdivided  to  include  the  budget  for  the  entire 
grid,  and  a  nested  budget  that  includes  only  the  microburst 
domain,  as  was  done  for  the  momentum  budget  studies.  The  eddy 
kinetic  energy  budget  may  provide  information  on  the  physical 
mechanisms  responsible  for  the  intensification  of  the 
raicroburst/misocyclone  system. 


6.5.2._1  The  Eddy  Kinetic  Energy  Budget  at  1845  MDT 

The  1845  MDT  case  (Table  10)  reveals  that  the  dominant 
terms  in  the  budget  are  the  HFD,  VFD ’  and  HFD '  terms.  The  total 
buoyancy  (BUP  +  PDG)  is  also  a  major  term  in  the  budget.  The 
generation  term  (WPG)  is  significant  at  levels  above  0.75  km, 
while  frictional  dissipation  (FDI)  is  a  major  contributor  at 
the  lowest  layer.  The  microburst/misocyclone  can  be  shown  as  a 
mechanism  for  transferring  energy  from  the  turbulent  state  to 
the  mean  state  for  this  particular  raicroburst  event.  The  inten¬ 
sification  of  the  microburst  might  seem  to  predict  a  movement 


of  energy  from  the  mean  to  the  eddy  state,  but  recall  that  the 


energy  source  is  believed  derived  from  the  misocyclone  level. 


The  mean  to  eddy  transformation  (term  WRS  >  0)  is  more  likely 


to  occur  with  the  misocyclone  than  with  the  surface  layer 


raicroburst.  At  levels  above  1.0  km,  there  is  a  increasingly 


positive  gain  in  eddy  kinetic  energy  through  the  VFD'  term,  as 


energy  is  extracted  from  the  mean  state  via  the  misocyclone.  A 


corresponding  loss  occurs  in  the  VFD  term,  representing  the 


mean  flow.  The  mean  kinetic  energy  is  extracted  from  the  high 


speed  flow  at  upper  levels  by  convergence  and  entrainment 


processes  within  the  misocyclone.  At  lower  levels  the  energy 


contained  within  the  eddy  flow  descending  within  the  microburst 


is  converted  back  into  the  mean  flow  at  low  levels.  Thus,  the 


microburst  tends  to  balance  the  difference  in  energy  levels 


between  the  top  and  bottom  of  the  misocyclone/microburst  layer. 


As  the  storms  intensify,  the  generation  term  serves  as  a  sink, 


caused  by  an  inability  of  the  mass  field  to  adjust  to  the  wind 


field.  The  high  speed  air  brought  down  from  aloft  enhances  the 


flow  from  the  low  pressure  at  upper  levels  to  the  high  pressure 


associated  with  the  surface  stagnation  high.  This  process  des¬ 


troys  kinetic  energy.  The  total  buoyancy  term  is  dominated  by 


the  precipitation  drag  component.  Precipitation  loading  is  an 


important  factor  in  the  maintenance  of  the  eddy  kinetic  energy 


budget,  and  in  maintaining  the  strong  downdraft.  The  friction 


term  is  important  in  the  lowest  layer,  where  it  is  a  relatively 


m 


large  sink  of  kinetic  energy.  It  remains  a  sink  above  the 


lowest  layer,  although  it  is  two  orders  of  magnitude  smaller. 
It  increases  at  upper  levels.  This  is  associated  with  the  tur¬ 
bulent  eddy  motion  and  shears  associated  with  the  roisocyclone 
at  those  levels.  This  result  is  consistent  with  the  observed 
maximum  of  the  turbulent  components  of  the  momentum  budget.  The 
dominant  mechanism  present  in  the  eddy  kinetic  energy  budget  is 
the  flux  divergence  of  kinetic  energy  from  the  mean  to  the  eddy 
flow  in  the  misocyclone,  and  from  the  eddy  to  the  mean  flow  in 
the  outflow  surface  layer  of  the  microburst.  The  intensifica¬ 
tion  of  the  misocyclone/microburst  circulation  at  1850  MDT  will 
be  examined  to  see  if  any  changes  occur  in  the  eddy  kinetic 
energy  budget. 


6 .5 .2^.  2  The  Eddy  Kinetic  Energy  Budget  at  J  350  MDT 
without  Tendency 

The  1850  MDT  case  (Table  11)  is  similar  to  the  1845  MDT 
case.  The  flux  divergence  is  still  the  predominate  mechanism 
with  precipitation  drag  also  important.  Several  important 
differences  are  noted.  The  first  of  these  is  the  increase  in 
the  magnitude  of  the  flux  divergence  terms  and  the  generation 
terra.  The  large  value  of  kinetic  energy  destruction  present  at 


1845  MDT  is  more  pronounced  at  1850  MDT.  The  lack  of  strong 
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horizontal  pressure  patterns  and  associated  eddy  circulations 


at  0.75  km  results  in  a  dominance  of  the  generation  term  by 
vertical  interaction  of  the  wind  field  with  the  pressure  field. 
Downward  flow  toward  high  pressure  at  0.75  km  will  thus  produce 
a  destruction  of  energy.  Enhanced  downward  flow  also  appears  to 
increase  the  magnitude  of  the  precipitation  drag  term.  Enhanced 
downward  motion  with  a  large  value  of  precipitation  drag  serves 
to  enhance  the  spin  up  of  the  misocyclone  through  vertical  vor- 
ticity  stretching.  It  also  suggest  a  lowering  of  the  reflec¬ 
tivity  cores  within  the  microbursts.  The  spin  up  and  compres¬ 
sion  of  the  circulation  allows  VFD'  to  achieve  a  positive  value 
even  at  0.50  km.  The  horizontal  flux  divergence  of  eddy  kinetic 
energy  achieves  a  large  value  of  dissipation  at  levels  below 
0.75  km.  The  increased  outflow  associated  with  the  intensified 
microburst  is  reflected  in  the  large  value  of  the  horizontal 
flux  divergence  of  eddy  kinetic  energy.  A  large  increase  in 
negative  buoyancy  at  1.25  km  reflects  the  rapid  intensification 
of  the  M2  misocyclone.  The  M2  misocyclone  is  warm-cored  at  1.25 
km,  since  the  mixing  of  dry  air  by  misocyclone  entrainment  is 
not  yet  sufficient  to  create  a  cold  core  downdraft.  The  sinking 
of  warm  air  creates  a  negative  value  of  buoyancy.  The  dominant 
physical  mechanism  associated  with  the  microburst  is  the 
extraction  of  mean  energy  by  the  misocyclone  and  transporting 
it  downward  to  the  raicroburst  level.  The  raicroburst  transfers 


eddy  kinetic  energy  back  to  the  mean  flow.  It  appears  that  the 
larger  the  loss  of  mean  flow  at  tne  misocyclone  level,  the 
greater  the  intensification  of  the  entire  circulation.  At  1850 
MDT  the  conversions  of  kinetic  energy  are  not  as  large.  Perhaps 
the  greatest  intensification  of  the  system  has  occurred  between 
1845  and  1850  MDT.  An  examination  of  a  separate  case  with  ten¬ 
dency  included  in  the  perturbation  pressure  calculation  is 
included  for  rough  comparison  purposes  only. 


6. 5. 2. 3  The  Eddy  Kinetic  Energy  Budget  at  1850  MDT  with 


Tendency 


The  inclusion  of  tendency  has  produced  profound  changes  in 
the  BUP  and  WPG  terms  of  the  equation  (Table  12).  The  genera¬ 
tion  term  retained  the  sign  of  the  previously  calculated  WPG 
term  for  the  non-tendency  case.  The  magnitude  of  the  term  has 
increased  by  a  factor  of  three  to  five.  The  buoyancy  term 
changed  sign  to  a  positive  value  and  increased  in  magnitude  by 
a  factor  of  three  to  five.  The  result  suggests,  for  a  rapidly 
evolving  complex  microburst  event,  the  pressure  calculation 
must  include  the  tendency  term.  The  inclusion  of  tendency  will 
affect  pressure  gradients,  which  in  turn  affects  those  terms  of 
the  equation  that  are  a  function  of  the  pressure  gradient.  The 
exact  nature  of  the  difference  and  its  physical  significance 


w, 


will  not  be  examined  in  this  study.  This  type  of  work  should  be 


accomplished  for  a  case  that  contains  three  or  more  time 
periods  at  two  minute  intervals  to  maximize  temporal  resolu¬ 
tion. 


6. 5. 2. 4  The  Eddy  Kinetic  Energy  Budget  at  1649  MPT 

The  final  case  examined  will  be  the  eddy  kinetic  energy 
budget  for  the  simple  case  microburst.  The  table  is  not  shown 
because  of  the  relative  simplicity  of  the  result.  The  domain 
considered  is  the  microburst  domain  (4X4  km).  The  balance  is 
almost  exclusively  between  the  flux  divergence  terms  and  the 
buoyancy  term.  The  most  interesting  difference  is  that  the  sim¬ 
ple  case  shows  a  conversion  of  mean  to  eddy  flow  at  all  levels 
and  of  a  similar  magnitude.  It  also  shows  a  gain  of  eddy 
kinetic  energy  at  the  lowest  levels  with  a  deficit  at  levels 
above  1.0  km,  which  is  a  similar  pattern  to  the  complex  case.  A 
common  characteristic  of  these  two  very  different  microburst  is 
the  role  that  eddy  kinetic  energy  flux  divergence  plays. 

The  vertical  flux  profiles  of  eddy  kinetic  energy  and  the 
eddy  kinetic  energy  budgets  for  1845  and  1850  MDT,  5  August 
1982  have  identified  the  raisocyclone  as  the  structural  mechan¬ 
ism  responsible  for  the  intensification  of  the  M2  microburst. 
The  eddy  kinetic  energy  budget  confirmed  the  key  role  played  by 


192 


^ww'wwwiwiwuii 


the  misocyclone  in  the  microburst  intensification  by  extracting 
energy  from  the  mean  flow  and  transporting  it  downward  through 
the  downdraft  to  low  levels,  where  the  horizontal  outflow  con¬ 
verts  the  eddy  kinetic  energy  transported  downward  from  aloft, 
to  the  low  level  mean  flow.  The  results  of  the  kinematic 
analysis,  the  momentum  budgets,  and  the  eddy  kinetic  energy 
budgets  may  be  integrated  with  results  of  previous  research  to 
provide  a  comprehensive  conceptual  model  of  the  mechanisms  and 
interactions  that  have  occurred  within  the  complex  and  simple 
raicroburst  domains. 


that  transports  the  air  further  downward.  The  stronger  environ¬ 
mental  wind  present  assures  that  the  environmental  scale  will 
dominate  the  raeso  y-scale  of  the  microburst.  No  well-organized 
raisocyclone  is  present,  because  of  the  dominance  of  the  larger 
scale  winds  and  the  lack  of  directional  shear  in  the  environ¬ 
mental  winds.  Srivastava  (1987)  has  shown  that  microbursts  may 
be  initiated  and  maintained  by  dry  entrainment  processes.  The 
key  to  this  type  of  microburst  seems  to  be  the  flux  rate  of 
high-moraentum,  dry,  environmental  air  across  the  concave  back 
edge  of  the  bow  echo.  The  greater  the  flux,  the  greater  the 
mixing  rate,  which  increases  the  downward  transport  of  air.  The 
rapidly  descending  air  is  compressionally  warmed  to  create  a 
warm  core  microburst,  because  of  vertical  motion  in  excess  of  5 
m/s.  This  is  a  sufficiently  rapid  descent  to  not  allow  the 
saturated  downdraft  to  mix  with  the  environmental  air.  The  air 
descends  nearly  dry  adiabatically  instead  of  moist  adiabati- 
cally.  These  results  were  produced  or  observed  in  a  simple 
one-dimensional  model,  Srivastava  (1985,1987);  PAM  station 
observations  Fuj ita  (1985);  and  in  a  study  using  dual-Doppler 
data,  Lin  and  Hughes  (1987).  The  microburst  centers  were  found 
to  be  centers  of  high  perturbation  pressure.  The  simple  case 
is  simple  in  part,  because  the  larger  scale  flow  can  dominate 
the  small  scale  flow;  therefore,  eddy  motions  of  misocyclone 
circulations  are  suppressed  or  are  of  secondary  importance,  the 


microburst  is  isolated  with  interactions  limited  to  the  outer 


edges  of  the  microburst,  and  these  types  of  storms  can  be 
predicted  by  techniques  based  on  the  presence  of  mid-level 
moisture  and  a  deep,  dry,  well-mixed  layer,  in  the  manner  of 
Caracena  et  al .  (1983).  Of  primary  interest  is  what  are  the 
driving  forces  behind  the  complex  case  of  a  multiple  microburst 
case  that  occurs  in  a  microburst-producing  storm  not  predicted 
by  present  forecasting  techniques. 

The  complex  case  is  characterized  by  several  microbursts, 
within  a  small  domain.  These  storms  interact  with  each  other  to 
produce  gust  fronts,  and  rotor  circulations,  and  they  distort 
each  other's  outflow  patterns.  There  are  similarities  with  the 
simple  case.  The  bow  echoes  do  develop  after  the  onset  of  the 
surface  raicroburst.  These  bow  echoes  are  aligned  perpendicular 
to  the  mean,  environmental,  high-momentum  air,  associated  with 
a  low  level  wind  maximum.  The  microbursts  are  coincident  with 
the  high  reflectivity  cores.  Inflow  into  the  bow  echo,  at  lev¬ 
els  above  1  kra,  is  mixed  with  the  saturated  downdraft.  Subse¬ 
quently,  the  high- momentum  wind  is  collected  above  1  km,  then 
carried  down  to  the  surface  layer  of  the  microburst.  The  role 
of  transport  and  buoyancy  processes  seems  common  to  both  types 
of  microbursts.  The  raicroburst  in  a  general  conceptual  model 
may  be  thought  of  as  a  mechanism  for  transporting  mass,  energy, 
and  momentum,  between  the  boundary  layer  (z  <  1  km)  and  the 


subc loud/parent  storm,  which  occurs  above  1  km  (AGL).  A  micro¬ 
burst  is  simply  a  meso  -scale,  organized,  turbulent  eddy.  The 
structural  features  responsible  for  the  microburst  development, 
and  the  physical  processes  that  maintain  them  are  not  unique. 

The  structural  features  for  the  complex  case  feature  a 
turbulent  upper  layer  above  1  km,  dominated  by  a  misocyclone 
circulation,  a  weakly  divergent,  weakly  eddying,  mid-level  at 
0.75  km  that  was  dominated  by  the  mean  flow,  and  a  very  tur¬ 
bulent  surface  layer  below  0.75  km,  which  was  dominated  by  the 
microbursts  and  their  secondary  circulations.  Of  particular 
interest  is  the  structure  and  development  of  the  misocyclone  at 
upper  levels,  and  the  microburst  circulation  at  low  levels. 

The  structure  of  the  storm  can  be  understood  by  examining 
how  the  microburst  may  have  developed.  The  complex  case  is  not 
conducive  to  forecasting,  using  techniques  based  on  negative 
buoyancy  from  a  deep,  dry,  mixed  layer  (Caracena  et  al . 
(1983)).  Dry  entrainment  and  evaporative  cooling  are  not  suffi¬ 
cient  to  initiate  the  raicrobursts  for  this  type  of  storm.  Pre¬ 
cipitation  loading  is  not  the  likely  cause,  because  tropical 
and  subtropical  convection  features  heavy  rains  with  compara¬ 
tively  mild  horizontal  surface  winds.  Hail  has  been  noted 
within  the  core  of  raicrobursts,  even  within  thunderstorms 
occurring  in  the  moist  gulf  coast  region.  Wakimoto  (1988)  used 


dual-polarization  data  from  the  CP-2  radar  to  locate  a  65  dBz 
hail  shaft  within  the  microburst  core.  He  suggested  that  the 
presence  of  a  hail  shaft  may  be  indi..»:ive  of  the  presence  of  a 
raicroburst.  A  study  of  a  Colorado  supercell  using  the  CP-3  and 
CP-4  radars  was  conducted  by  Miller  et_  al.  (1988).  The  results 
indicated  that  the  location  of  hail  reports  were  mostly  inside 
the  swath  of  the  55  dBz  contour.  The  hail  sizes  noted  were  2-3 
cm  in  diameter  with  a  maximum  reflectivity  of  65-75  dBz.  Using 
modeled  data,  Foote  (1984)  showed  that  hailstones  of  1.0  cm  in 
diameter  lose  73Z  of  their  mass  during  descent  through  warm 
air.  Many  researchers,  including  Hjelmfelt  (1984)  and  Parsons 
et  al.  (1985)  have  noted  that  cyclonic  circulations  aloft  are 
associated  with  appendages  of  very  high  reflectivity.  Parsons 
et  al .  (1885)  indicated  that  for  a  30  June  1982,  JAWS  case, 
the  raicroburst  downdraft  initially  formed  close  to  a  small  core 
of  reflectivity  embedded  within  a  larger  'arm'  of  reflectivity. 
An  intense  cyclonic  circulation  centered  about  a  reflectivity 
appendage,  with  continuity  to  heights  above  cloud  base,  was 
noted  by  Hjelmfelt  (1984).  He  also  noted  that  peak  surface 
reflectivities  were  attained  prior  to  the  maximum  microburst 
divergence  intensity.  Also,  evidence  of  low  relectivity  notches 
in  the  neighborhood  of  the  raicroburst  downdrafts  suggested  that 
melting  provides  an  additional  source  of  energy.  From  a  three- 
dimensional  model,  Proctor  (1988)  showed  that  a  warm  core  may 


be  present.  Hail  is  too  small  to  compensate  for  compress ional 
heating,  unless  the  vertical  velocity  of  the  downdraft  is  very 
small.  Since  hail  plays  a  relatively  minor  role  in  maintaining 
the  microburst,  its  role  may  be  in  a  different  capacity.  The 
main  problem  in  distilling  the  answer  from  the  evidence  is  to 
determine  if  the  presence  of  the  raicroburst  is  the  result  of 
the  hail  or  the  hail  is  the  result  of  the  presence  of  a  raicro¬ 
burst.  Hail  falling  in  strong  downdrafts  has  a  significant 
advantage  over  hail  that  does  not,  in  that  it  can  spend  consid¬ 
erably  less  time  in  the  warm  air  before  reaching  the  ground.  A 
comparison  of  the  1649  MDT,  14  July  1982  and  the  1845  and  1850 
MDT,  5  August  1982  cases  shows  that  the  cloud  base  is  higher 
for  the  former  case.  The  lower  reflectivities  encountered  in 
the  simple  case  suggests  that  the  hail  was  incidental  to  the 
microburst.  However,  the  complex  case  occurred  in  an  environ¬ 
ment  of  higher  available  latent  heat  for  melting.  Also,  the 
downward  motion  within  the  microburst  was  one-half  of  the  mag¬ 
nitude  experienced  in  the  simple  case.  This  suggests  that  hail 
may  be  responsible  in  part,  for  the  microburst,  since  a  larger 
volume  of  hail  is  necessary  to  maintain  the  high  reflectivity 
values  in  the  complex  case,  where  the  melting  rate  is  higher, 
due  to  a  slower  descent  of  the  hail  in  the  microburst  down- 
draft,  relative  to  the  simple  case,  and  the  greater  availabil¬ 
ity  of  environmental  latent  heat.  If  there  was  another  causal 


mechanism  present,  the  reflectivity  core  should  have  been  lower 
or  the  same  as  the  simple  case.  I  believe  that  a  main  causal 
mechanism  for  the  complex  case  microburst  was  hail  melt 
processes,  but  dry  entrainment  processes  associated  with  the 
misocyclone  circulation  and  downward  momentum  transport  of 
energy  extracted  from  the  mean  state  to  the  eddy  state  main¬ 
tained  the  microburst  at.  lower  levels.  New  forecasting  tech¬ 
niques,  to  consider  the  role  played  by  hail,  are  necessary  to 
forecast  raicrobursts  for  relatively  moist  soundings.  By  examin¬ 
ing  the  conceptual  model  in  a  temporal  sense,  the  microburst 
synthesis  is  clearer. 

It  is  postulated  that  the  evaporative  cooling  from  precip¬ 
itation  initiates  the  downdraft,  see  Srivastava  (1985),  but  the 
focus  for  the  raicroburst  is  provided  by  the  descending  hail 
shaft.  The  hail  melt  contributes  the  necessary  additional  nega¬ 
tive  buoyancy,  when  combined  with  evaporative  cooling,  to 
create  a  sufficient  downward  flux  of  mass  and  momentum  to  pro¬ 
duce  the  microburst.  The  hail  melts  below  the  freezing  level, 
with  the  maximum  contribution  of  hail  melt  to  negative  buoyancy 
occurring  at  the  bottom  of  the  descending  microburst.  At  that 
level,  the  melting  is  maximized,  as  the  cool  microburst  air 
mixes  with  the  warm  subcloud  layer.  A  commonly  observed  event 
with  raicrobursts  is  a  descending  reflectivity  core,  accompanied 
by  a  'melting  level  depression'.  The  contribution  from  hail 
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melt  to  the  microburst  diminishes,  once  the  microburst  impinges 
upon  the  surface.  The  air  through  the  subcloud  layer  has  been 
cooled;  therefore,  melting  slows  down  with  time.  This  is  con¬ 
firmed  by  the  descending  reflectivity  core  and  the  melting 
level  depression.  As  the  hail  shaft  becomes  more  insulated,  the 
hail  melt  contribution  to  negative  buoyancy  is  limited  to  the 
outer  edge  of  the  reflectivity  core.  Since  hail  melt  no  longer 
becomes  a  major  factor  in  the  microburst's  development,  some 
other  mechanism  must  be  evolving  to  take  its  place. 

The  misocyclone  is  not  present  during  the  early  stages,  as 
confirmed  by  M3  in  the  5  August  1982  case.  As  the  downdraft  and 
microburst  progress  downward  toward  the  surface,  the  atmos¬ 
pheric  column,  through  which  the  downdraft  descends,  becomes 
stretched.  The  stretching  occurs  in  a  light  wind  environment  of 
environmental  winds  less  than  10  m/s,  but  with  strong  direc¬ 
tional  shear  of  the  winds  backing  with  height.  Since  the  raeso 
/-scale  dominates  the  raeso  -scale  winds,  the  stretching  of  the 
vorticity  initiates  a  cyclonic  gyre  above  0.75  km.  The  gyre 
entrains  dry,  high-raoraentum,  environmental  air  into  the  down- 
draft,  where  it  mixes  with  the  cool  saturated  downdraft.  The 
enhanced  dry  entrainment  with  evaporative  cooling,  downward 
flux  of  momentum,  and  enhanced  precipitation  drag  due  to  the 
precipitation  core,  become  the  primary  raaintainers  of  the 
microburst.  The  weak  reflectivity  notches,  noted  in  this 


study,  are  a  manifestation  of  this  process.  The  evolution  of 
the  bow  echo  parallels  the  development  of  the  misocyclone.  In 
this  study  the  rotating  head  of  the  bow  echo  (Fujita  (1985)) 
was  coincident  with  the  steady  state  of  the  1649  MDT  case,  and 
Ml  of  the  5  August  1982  case. 

The  surface  raicroburst  is  a  direct  result  of  the  downdraft 
impinging  upon  a  horizontal  surface.  The  stretched  downdraft, 
becomes  compressed  as  its  downward  velocity  is  slowed  from  con¬ 
tact  with  the  ground.  The  compression  results  in  a  spin  down  of 
the  downdraft.  The  cyclonic  circulation  undergoes  spin  down  to 
a  neutral  condition  at  0.75  km.  Below  that  level,  a  rapid  spin 
down  creates  a  highly  divergent  circulation.  The  microburst 
sets  up  secondary  circulations  on  its  margins.  A  cat’s  eye  roll 
vortex  forms  on  the  outer  edge  of  the  microburst,  when  the  hor¬ 
izontal  shear  between  the  raicroburst  outflow  opposes  each 
other.  The  higher  density  flow  will  determine  the  sense  of 
rotation.  Interieri  et  al_.  (1988)  used  Doppler  lidars  to  show 
that  uneven  flows  produce  gyres  and  cat's  eye  circulations.  In 
addition  the  gust  front  and  environmental  flow  combine  to 
create  an  updraft  that  stretches  the  surface  vorticity  to 
create  a  spinning  updraft.  At  misocyclone  level,  the  misocy¬ 
clone  and  mesocyclone-like  circulation  interact.  The 
mesocyclone-1 ike  circulation  facilitates  the  entrainment  of 
environmental  air  into  the  misocyclone,  thus  becoming  a 
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feedback  mechanism  within  the  misocyc lone/microburst  circula- 
t  ion . 


Finally,  a  summary  of  the  main  points  of  the  conceptual 
model  and  the  microburst  structure  are  presented  as  follows: 
imicroburst  structure  are  presented  as  follows:  1)  Fujita's 
(1985)  criteria  of  the  microburst  are  to  be  updated.  The  new 
criteria  should  include  the  downward  vertical  velocity,  reflec¬ 
tivity,  and  the  maximum  horizontal  wind  shear;  2)  the  role  of 
the  misocyclone  circulation,  occurring  above  0.70  km,  is  to  be 
investigated;  7)  the  role  of  a  circulation  occurring  where  the 
environmental  wind  and  the  gus'  front  wind  shear  are  a  maximum, 
should  be  further  researched.  *s  other  microbursts  affect  the 
infraction  between  storm  and  environment,  the  position  of  the 
mesocyclone-like  cyclonic  circulation  shifts  to  the  new  area  of 
highest  shear.  4)  The  mesocylone-like  circulation  interacts 
with  the  misocyclone  to  feed  dry  environmental  air  into  the 
misocyclone.  5)  The  raicroburst  occurs  before  the  creation  of 
the  misocyclone  and  the  bow  echo.  The  misocyclone  and  the  bow 
echo  intensify  in  unison.  The  rotation  head  of  the  echo  is 
coincident  with  the  misocyclone.  6)  The  misocyclone  assists  in 
entraining  and  mixing  the  dry  environmental  air  with  the 
saturated  downdraft.  7)  The  inflow  into  the  downdraft/bow  echo 
occurs  at  the  concave  back  edge.  This  edge  is  aligned  perpen¬ 
dicular  to  the  environmental  flow.  8)  The  downdraft,  for  the 


wet  Colorado  microburst,  appears  to  originate  from  a  reflec¬ 
tivity  core  of  at  least  45  dBz.  9)  The  momentum  flux  and  eddy 
kinetic  energy  flux  were  downward  within  the  raicroburst.  This 
indicates  a  net  transport  of  mass  and  momentum  to  the  surface 
from  aloft.  10)  Frictional  parameterization  and  the  inclusion 
of  the  local  tendency  term  improves  the  momentum  check  (Table 
13).  11)  The  microburst  center  is  an  area  of  high  perturbation 
pressure,  with  an  area  of  relative  low  pressure  surrounding  the 
microburst.  12)  Roll  vortices  occur  at  the  microburst's  edge.  A 
maximum  of  vortice  generation  occurs  between  colliding  outflows 
from  separate  microbursts,  with  a  secondary  maximum  where  the 
environmental  wind  and  the  microburst  outflow  intersect  at  a 
wind  discontinuity  of  180  deg..  13)  For  downward  motion  of  less 
than  3-4  m/s,  a  cold  core  microburst  results.  If  the  downdraft 
is  of  a  greater  magnitude,  a  warm  core  will  result. 

The  conceptual  model,  derived  from  this  study,  is  based  on 
a  composite  of  the  present  research  and  findings  from  previous 
work.  Additional  work  is  required  to  confirm  or  modify  this 
model,  and  to  answer  questions  implied  by  the  present  research 


but  not  answered. 
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8.  Suggested  Topics  for  Future  Research 

Much  work  needs  to  be  conducted  to  verify  or  modify  these 
results.  The  following  is  a  list  of  possible  topics  for  future 
study:  1)  A  trajectory  study  using  data  from  this  case  and 
similar  cases.  2)  A  study  of  pressure  retrieval  using  a  micro¬ 
physical  parameterization  of  hail  melt  should  be  conducted  to 
determine  if  it  may  or  may  not  have  a  significant  impact  on  the 
perturbation  pressure  recovery  and  the  perturbation  temperature 
recovery.  3)  A  study  using  dual-Doppler  data  could  be  con¬ 
ducted,  using  a  two  minute  interval,  centered  difference,  velo¬ 
city  tendency  term  for  calculating  the  perturbation  pressure. 
4)  A  numerical  model  should  be  used  to  confirm  or  modify  the 
conceptual  model.  5)  A  study  using  a  0.25  km  grid  spacing 
should  be  conducted  to  examine  the  impact  of  the  inherent 
smoothing  of  the  Barnes's  (1973)  scheme  on  the  small  scale 
features.  6)  A  study  using  a  0.125  km  grid  level  (AGL)  may 
show  more  detailed  structure  of  the  storm  at  the  lowest  level, 
and  a  more  accurate  result  from  the  upward  integration  of  the 
anelastic  continuity  equation.  7)  Several  data  sets  from  dif¬ 
ferent  geographical  regions  should  be  used  to  conduct  similar 
studies.  8)  A  study  on  how  sensitive  the  individual  components 
of  the  budget  equations  are  to  random  errors  should  be  con¬ 
ducted.  9)  Additional  work  on  the  5  August  1982  case  should 
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include  the  1847  MDT  time  period  to  allow  the  researcher  more 


insight  into  the  microburst’s  evolution.  10)  A  vorticity 
budget  should  be  conducted  to  see  how  stretching  of  the  vorti¬ 
city  influences  the  growth  and  decay  of  the  raisocyclone  and  the 
microburst.il)  A  forecasting  limiting  index  to  indicate  if  hail 
is  necessary  to  generate  a  microburst,  and  to  what  degree  would 
that  be  indicated  by  examining  the  reflectivity  return,  should 
be  developed.  12)  A  correlation  study  between  observed  reflec¬ 
tivity  values,  the  wet  bulb  of  the  subcloud  layer,  and  the 
raicroburst  intensity  should  be  conducted,  as  a  method  for  fore- 


Appendix  A:  Frictional  Parameterization 


r 


Friction  nay  be  conceptually  thought  of  as  cellular  mixing 
circulations  that  efficiently  transmit  energy  between  areas  of 
deficit  and  surplus,  principally  down-gradient.  Except  at  the 
microscale,  below  several  millimeters,  where  molecular  dissipa¬ 
tion  is  significant,  mesoscale  and  synoptic  scales  are  dom¬ 
inated  by  eddy  transport  of  momentum,  heat  and  moisture.  Since 
eddy  circulations  are  two-dimensional  eddies,  the  degree  of 
correlation  between  the  magnitudes  of  perturbations  of  the 
three-dimensional  components  of  the  mean  wind  indicates  the 
strength  and  efficiency  of  the  frictional  eddies.  The  eddy 
strength  is  stated  symbolically  and  may  be  calculated  directly, 
if  data  response  time  is  sufficiently  rapid,  as: 


where  p  is  the  mean  air  density,  and  i  and  j  vary  from  1  to  3. 
This  method  produces  a  qualitative  field  of  eddy  stresses, 
which  describes  the  momentum  flux  exchange  conducted  by  the 
eddies,  but  are  sensitive  to  error  for  quantitative  calcula- 
t ions . 

The  quantitative  method  used  to  calculate  the  ^ij  and  Km 
values  for  the  eddy  momentum  and  kinetic  energy  budget  begins 


Eq.  (2) 


with  the  Reynold's  equations  of  motion  in  tensor  form: 

where  the  overbar  represents  the  mean  of  the  quantity,  the 
prime  represents  an  anomaly,  eddy  motion,  £ijk  is  a  unit  ten- 
sor.  L  ij  is  the  molecular  stress  tensor  and  other  variables 
are  by  standard  convention.  Substituting  Eq.  (1)  into  Eq.  (2) 


(3) 


results  in 

\ZEj.  -  '£■  CPU.U: 

Jt  .  ' 

A  B  C  D 

where  Term  A  «  local  time  rate  of  change  in  mean  momentum,  per 


unit  volume,  Term  B  »  flux  convergence/divergence  of  momentum 


by  mean  motion,  Term  C  -  Coriolis  force,  per  unit  volume.  Term 
D  -  geopotential  gradient  force,  per  unit  volume,  and  Term  E  « 
flux  convergence/divergence  of  momentum  by  eddy  and  molecular 
motion. 


The  pressure  gradient  force  is  contained  within  term  E  due 
to  the  molecular  stress  (Tij).  Klemp  and  Wilhelmson  (1978) 
parameterized  the  eddy  stress  gradient  of  the  subgrid  turbulent 
momentum  terms,  f.,  as: 


The  local  momentum  flux  may  be  defined  as  the  loss  of 


momentum  to  other  scales  of  motion  from  the  total  subgrid-scale 
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eddy  kinetic  energy  as  follows: 

Ci7u;  =  -KM(^.^;V|iijE  ei.  (5> 

where  Km  is  the  eddy  momentum  transfer  coefficient, 

ij  is  the  Kronecker  constant,  and  E  is  the  subgrid-scale  eddy 
kinetic  energy,  per  unit  mass,  symbolically  stated  as: 

E  =  i(^2)  Eq.  (6) 


Ideally  one  would  like  to  be  able  to  calculate  friction 
directly  with  fast  response  sensors  (Doppler  radar  is  not). 
Doppler  radar  data  are  actually  sample  volumes,  not  points, 
which  are  linearly,  temporally  interpolated.  Since  this  isn't 
conducive  to  direct  calculation  of  turbulance  at  the  subgrid 
scale  and  the  energy,  E,  contained  within  the  eddies,  the  fol¬ 
lowing  technique  is  used. 


To  solve  this  problem,  parameterization  is  used.  From 
Sraagorinsky  (1963)  we  get  a  form  of  Km  originally  used  to  cal¬ 
culate  Km  at  a  larger  scale  of  motion  from  the  deformation  ten¬ 
sor.  Deardorff  (1972)  proposed  that  for  unstable  environments, 
the  subgrid-scale  eddy  coefficient,  c,  is  0.21.  The  value,  c, 
is  an  efficiency  type  quantity  that  parameterizes  the  energy 
cascade  to  lower  scales.  The  result  is  as  follows: 


**  "  Tz  **il 


Eq.  (7) 

Knowing  Km  one  can  obtain  E  by  dimensional  grounds  as  pro- 


211 


and 


(On=  #X|2$1 


Eq.  (13) 


The  Reynold's  stresses  may  be  expanded  in  similar  fashion: 

rly -  T,’  ■■  ■-  7f%l£ 

Tu  - -  r {( (& .£)-fJ,,e}=  'V^  I S s)(. 

r„  •-  =  ?f  ^la*gl(g*S)J, 

r«.H*  «w>»ifc**0*i‘*ei 

-  "  £  1  it  »y  K  **  ^  / 


Eq.  (14) 


iil 


■*  1  i  \ .  iw  itr 


The  Reynold's  stresses  were  calculated  from  Eq.  (15)  and  were 
utilized  in  the  pressure  calculation  and  the  turbulent  kinetic 


energy  budget. 


Appendix  B:  Additional  Error  Analysis 

The  following  scale  analysis  was  calculated  to  show  the 
impact  of  possible  error  and  uncertainty  of  hail  on  the  three- 
dimensional  wind  field.  This  process  is  important  not  only  to 
determine  the  validity  and  relevance  of  the  derived  fields,  but 
also  to  demonstrate  their  sensitivity  to  error. 

The  first  possible  source  of  error  that  will  be  examined 
is  the  impact  of  hail  and  attendant  reflectivity  uncertainty  on 
the  computation  of  the  three-dimensional  wind  field,  Armijo 
(1969).  The  complete  set  of  equations  was  derived  in  Section  3. 

Assume  a  point  very  close  to  both  radars  with  hail  occur¬ 
ring  with  a  terminal  fall  speed  of  20  ra/s.  In  a  shallow  volume 

such  as  the  one  used  in  this  study,  the  elevation  angle  was 

o 

never  more  than  10  and  the  closest  point  to  a  radar  was  6.5 
km.  The  highest  level  analyzed  was  1.25  km  for  this  case,  the  5 
August  1982  case,  see  Fig.  (Bl).  The  result  from  calculation. 

tan  -  y/x  -  tan  10  deg.  -  y/20  m/s 
y  -  3.5  m/s 

This  condition  requires  the  entire  pulse  volume  to  behave 
as  one  hailstone.  Realistically,  a  size  distribution  for  hail, 
given  by  Auer  (1972),  gives  a  more  physically  meaningful  result 


\V- 


by  specifying  the  amount  of  hailstones  per  volume  of  air. 


$ 

ft 


C 


-34 

N(d)  -  561.3d 


Where  N  (mm  )  and  d  (mm)  for  1  <  d  <  50  mm.  A  stone  of  20 


m/s  fall  speed  would  correspond  to  a  stone  at  the  upper  end  of 


the  d  (mm)  scale.  As  the  relationship  suggests,  the  greater  the 


individual  hail  size,  the  greater  the  possible  contribution  to 


the  radial  velocity,  but  the  less  the  actual  number  of  such 


hail  per  pulse  volume.  The  result  would  suggest  that  for  a 


volume  of  7.07  X  10  m  ,  the  ensemble  contribution  would  be 


negligible  in  the  boundary  layer  unless  an  intense  hailstorm  is 


present . 


Of  greater  importance  is  the  impact  of  hail  on  the  termi-. 


nal  fall  speed  via  the  precipitation-drag  term.  The  relation¬ 


ship  chosen  was  the  Martner  (1975)  parameterization.  It  is 


weakly  dependent  on  reflectivity.  By  introducing  a  20X  error, 


which  seems  to  be  an  upper  limit  according  to  Spahn  and  Smith 


(1976),  results  as  follows: 


Error  V  (20X)  ~  .5  ra/s 


The  observed  reflectivity  was  chosen  to  be  56  dBz  with  an 


actual  value  of  70  dBz. 
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A  second  source  of  error  considered  here  is  the  impact 
error  in  terminal  fall  speed  may  have  on  the  horizontal  wind 
components.  Again,  an  extreme  example  is  used.  If  one  radar  is 
6.5  km  from  the  datum,  the  other  is  approximately  15  km  away, 
see  Fig.  (3.1).  The  radial  winds  are  an  average  of  8  m/s. 
Although  both  radars  will  not  likely  see  the  same  magnitude  of 
radial  velocity  or  even  the  same  sign  due  to  geometric  con¬ 
siderations,  they  will  for  purposes  of  illustration.  Plugging 
back  into  one  of  the  Armijo  (1969)  equations  and  using  the  same 
input  parameters  results  in: 


"flu 

tm 


VO 


'] 


u  -  tamo1  mV-DYi  -I(*q  nV-Oya- xipV> 

(V*,-x,ya) 
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,  or  substituting: 


It  is  obvious  from  this  illustration  that  the  precipita¬ 
tion  drag  term  will  introduce  a  less  than  10%  error  into  the 
calculation  even  when  assuming  the  presence  of  hail  and  a  max¬ 
imum  error  in  the  reflectivity.  Also,  the  Barnes  (1973)  routine 
will  smooth  out  any  areas  of  maximum  error  so  the  singular, 
exceptionally  large  values  of  error  are  minimized. 


The  same  procedure  with  more  typical  values  (reflectivity 


s 
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